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Abstract 

The publicly available spectrum generators for the NMSSM often lead to different predictions for the 
mass of the standard model-like Higgs boson, even though they all implement two-loop computations in 
the DR renormalization scheme. Depending on the parameter point, the differences can exceed 5 GeV, 
and even reach 8 GeV for moderate superparticle masses of up to 2 TeV. It is shown here that these 
differences can be traced back to the calculation of the running standard model parameters entering 
all calculations, to the approximations used in the two-loop corrections included in the different codes, 
and to different choices for the renormalization conditions and scales. In particular, the importance of 
the calculation of the top Yukawa coupling is pointed out. 


1. Introduction 

Now that run II of the Large Hadron Collider (LHC) has started, either the first clear hints for 
new physics beyond the standard model (BSM) will finally show up or many well-motivated theories 
for BSM physics will come under even greater pressure. While in run I a scalar was discovered 
mm which has all properties of the long-expected Higgs boson of the standard model (SM), at the 
same time stringent limits were placed on the simplest scenarios for new physics. In particular, the 
minimal realization of supersymmetry - the minimal supersymmetric standard model (MSSM) - has 
lost at least some of its appeal: no hints for supersymmetric (SUSY) particles were found, and the 
observed value of the SM-like Higgs mass requires either multi-TeV stop masses or special values of 
the stop-mixing parameter. However, if the stops are heavy and the MSSM is assumed to be a valid 
description of nature far above the electroweak scale, then a large fine-tuning is required to stabilize the 
vacuum expectation value (VEV) of the Higgs field. On the other hand, the large stop mixing needed 
to generate a Higgs mass of 125 GeV for stop masses around one TeV turns out to be dangerous 
because of charge and color breaking minima ED H ED E 0- This leads to the question of whether 
the MSSM is really a natural completion of the SM and has attracted more interest in other SUSY 
scenarios. This is especially the case for SUSY models that predict a larger Higgs mass already at tree 
level. The best-studied examples are models with singlet extensions, where the interaction between 
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the Higgs fields and the gauge singlet can raise the Higgs mass mmm and significantly reduce the 
fine-tuning 

mimnaiaisiiaimimimin]. AZs symmetry is often introduced to forbid all 
dimensionful parameters in the superpotential and to solve the fj, problem of the MSSM at the same 
time [9]. This realization of a singlet extension of the MSSM is the next-to-minimal supersymmetric 
standard model (NMSSM). The NMSSM comes with a rich collider phenomenology [2J] [921 [53], and 
also provides possibilities to hide SUSY at the LHC [24j . 

To confront the NMSSM with the measured Higgs mass, a precise calculation is necessary. The first 
calculations of radiative corrections to the Higgs masses in the NMSSM assumed an effective field theory 
below the SUSY scale (the scale of the stop/sbottom masses), and integrated the renormalization group 
equations (RGEs) for the corresponding couplings from the SUSY scale to the weak scale [251 [551 27 . 
One-loop corrections ^ of 0(at + otb) in the effective potential approximation have been computed in 
Refs. [531 221 GUI El ■ Two-loop corrections of 0(a s at + af) have been included in the NMSSM in 
the leading-logarithmic approximation (LLA) in Refs. [321133], and one-loop corrections to the SM-like 
Higgs mass including NMSSM-specific and electroweak couplings in the LLA in Ref. [34- Full one- 
loop calculations including the momentum dependence have been performed in the DR renormalization 
schem^] [35] [37] , or in a mixed on-shell (OS)-DR scheme [38} 139]; at the two-loop level, the dominant 
G(a s (ab + ctt)) corrections in the DR scheme have become available a few years ago [35]■ The two-loop 
corrections involving only superpotential couplings such as Yukawa and singlet interactions WL and 
a two-loop calculation of the 0(a s at) corrections with the top/stop sector renormalized either in the 
OS scheme or in the DR scheme mi, have also been provided recently. 

At the moment, five public spectrum generators for the NMSSM are available which make use 
of these results to compute the Higgs-boson masses and mixing matrices in the NMSSM for a given 
parameter point: the stand-alone codes NMSSMTools 0213311113, SOFTSUSY [45, 35, 37, and NMSSMCALC 
[48], as well as NMSSM versions of FlexibleSUSY [49] and SPheno [501 [51] . which are based on SARAH 
[55115311531 [551 [55] , The various tools, however, often make predictions that differ by several GeV for the 
mass of the SM-like Higgs boson even if using the same renormalization scheme. These differences are 
bigger than the ones observed between codes that calculate the Higgs mass in the MSSM [57][53]. The 
aim of this work is to identify the origin of these deviations for the DR calculation, and to determine 
the numerical significance of each difference between the codes. A comparison between the DR and 
OS calculations, which would provide an estimate of the theoretical uncertainty, is beyond the scope 
of this analysis and will be given elsewhere. 

The paper is organized as follows. In sec. [2] we introduce the model and the conventions used. In 
sec. [3] the different codes and their features are presented, and the results obtained in six test points 
when using the codes out-of-the box are shown. The differences in the Higgs-mass calculation between 
these codes are discussed in sec. [4] Here, we concentrate on a SUSY scale input, while we address 
additional effects in a GUT scenario in sec. [5] We conclude in sec.[Gl 


2. The NMSSM and its conventions 

We consider in the following a model with the particle content of the MSSM extended by a gauge 
singlet superfield S. The general, renormalizable and R-parity conserving superpotential for this model 
in terms of the (hatted) doublet ( H u ,H d ) and singlet superfields reads 

Wg = WVukawa + 4* (/i + XS)H u H d + —f-l s S 1 2 + t s S , (1) 

with 

WVukawa = YuQUHu + Y d QDH d + Y e LEH d , (2) 


1 We define c*t,b,T = ^tb t /(4tt), where Yt,b,r are the third-family Yukawa couplings, and adopt analogous definitions 
for a\ and where A and k, are the NMSSM-specific couplings in eq. 0 below. We follow the standard conven¬ 
tion of denoting, e.g., by (D(at) the one-loop corrections to the Higgs masses that are in fact proportional to 
i.e. proportional to Y t 4 . Similar abuses of notation affect the other one- and two-loop corrections. 

2 Strictly speaking, all calculations have been done in the DR^-scheme where the e scalars are decoupled m- 
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where we have suppressed color and isospin indices in every term. In general, the couplings in eq. ([I]) 
are complex parameters, and those in eq. ([2]) are complex 3x3 matrices. The most general soft 
SUSY-breaking terms for this theory are given by 

Koft = rn 2 s \S\ 2 + m 2 H jH u \ 2 + m 2 Hd \H d \ 2 + pmjf 

+ .MiAjj + M 2 X^ V + M 3 A“ + T d qdH d + T u quH u + T e elH d 

T K S 3 + h.c)j , (3) 

with the left-handed squark and slepton doublets / = q and l and the right-handed squarks and 
sleptons / = Ur, dn, Hr. All generation indices are suppressed. The soft SUSY-breaking masses m 2 , 
m 2 Hd and m 2 H ^ are always real, and the sfermion masses m 2 ~ are hermitian 3x3 matrices, while all other 
soft terms can be complex. The soft SUSY-breaking trilinears T\ and T K are often expressed in an 
alternative way, as the product of the corresponding superpotential coupling times a free parameter 
with the dimensions of a mass 


1 


1 


+ B,H u H d + -B s SS + Xs S + T\SH u H d + - 


T\ = AA> , T k = kA | 


(4) 


Similarly, the trilinear soft SUSY-breaking interactions involving sfermions can be expressed as 

Ti =YiAi , i = u,d,e . (5) 

It is common to introduce a Z 3 symmetry to forbid all dimensionful parameters (i.e., /r, fi s and t s ) 
in eq. Q. This solves the fj, problem of the MSSM m - at the price of introducing a domain wall 
problem j60[ - and leads to the NMSSM, where the surviving terms in the superpotential are 

WnMSSM = WVukawa + + A SH u H d , (6) 

and the soft SUSY-breaking terms BB s and \s in eq. ([3]) also vanish. 

For suitable values of the remaining s-dependent SUSY-breaking terms m s and T re , a VEV v s for 
the scalar singlet S of the order of the SUSY-breaking scale is generated. This dynamically generates 
the /i term which gives masses to the higgsinos, 


deS = MS) . (7) 

Note, there are two common conventions: (S) = v s or (S) = v s jy/2. The SUSY Les Houches Accord 
(SLHA) (EM does not fix a convention, and the codes compared here make different choices. However 
in this text we will use the latter convention. Finally, after electroweak symmetry breaking (EWSB) 
where also the neutral components of H d and H u receive VEVs v d and v u , the singlet mixes with 
the other doublets^ In the CP-conserving case, this leads to three CP-even physical scalars and two 
CP-odd ones together with the neutral would-be Goldstone boson. It is common to define the basis 
for the CP-even states as {(j) d , <f> u , <f> s ) where (j)j are the real parts of the neutral components of H d , H u 
and S. This leads to the tree-level mass matrix (M^) , derived from the tree-level scalar potential 


{Ml) 


tree 

*3 


^2^tree 

dfiidcfrj 


( 8 ) 


The pole masses (i = 1, 2,3), identified with the masses of the three CP-even physical scalars, are 

the eigenvalues of the momentum dependent, loop-corrected mass matrix 


A4V) = (A4 2 ) tree + ^nW( P 2 ), (9) 


3 While in the CP violating NMSSM the VEVs can be taken complex, it turns out that, after exploiting the minimum 
conditions of the scalar potential, the NMSSM Higgs sector at tree level can be described by one complex phase, which 
is given by (f y = <p K — ip\ + 2ip a — where ip u and (p s are the complex phases of the VEVs of H u and S, respectively, 
and and ip\ the complex phases of k, and A. 
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each evaluated at p 2 = m \.. The II^”^ ( p 2 ) denote 3x3 matrices of ?r-loop self-energies calculated for an 
external momentum p 2 . As mentioned in the introduction, only the one-loop contribution, n ( 1 )(p 2 ), 
is completely known so far. For 11^ only partial results in the limit p 2 = 0 exist in the NMSSM. 
Different choices of the renormalization scheme lead to different results for the physical scalar masses 
if the perturbative series in eq. © is truncated at a given order n. The differences between the 
results are of 0 (n(" +1 )), and provide an estimate of the theoretical uncertainty of the calculation due 
to uncomputed higher-order corrections. The matrix that diagonalizes M 2 {p 4 2 * ) is called Z H in the 
following© 

m d ™ 9 = Z H M 2 (jp 2 )Z T H . (10) 

In general, Zpj is a p 2 -dependent quantity. The NMSSM spectrum generators adopt different definitions 
for the loop-corrected rotation matrix Zh that they give as output: the external momentum is either 
set equal to zero, or set equal to one of the mass eigenvalues. 


We conclude this section by presenting the set of input parameters that enter the calculation of 
the Higgs-boson masses and mixing matrices in the CP-conserving NMSSM. First, there are the three 
tadpole conditions that ensure the minimum of the Higgs potential at non-vanishing VEVs v u ,Vd,v s , 


dV eS 

d<t>i 


= 0 , 

min 


i = d, u, s , 


( 11 ) 


where V cH represents the effective scalar potential computed at the same perturbative order as the 
mass matrix in eq. (19]) . These conditions fix three parameters. The simplest choice is to solve them 
for m 2 , ttVjj and rnjj , and this is done in all the comparisons for every code we use. Since one 
combination of VEVs (namely, v 2 = v 2 + v 2 ) is constrained by experimentally measured observables, 
this leaves six free parameters in the Higgs sector of the model: 


A, k, A x , A k , tan/3, p eB , ( 12 ) 

where tan/3 = v u /vd, and p eB stands in for v s . Under the assumption that the soft SUSY-breaking 
terms are flavor-diagonal and that only third-generation soft interactions are non-negligible, this re¬ 
duces the additional input parameters required to fully define the model to the set 


A t , A bl A t , mi L , mp R , M 1 , M 2 , M 3 , Mj , ( 13 ) 

with m 2 ^ = m~ 33 , m = rn| 3 ^J A T = A 6 j 33 , A b = ^^ 33 , A t = A Uj 33 , and where we assumed for 
simplicity that all other diagonal entries of m 2 are degenerate and equal to Mj . 


One can further reduce the number of input parameters by assuming a specific mechanism to 
transmit SUSY breaking from a hidden sector to our visible sector. In models inspired by minimal 
supergravity the following relations among the soft terms exist at the gauge-coupling unification scale: 

m 2 ~ = mil, Mi = M 1/2 (i = 1, 2 , 3) , Aj = A a ( j = e,d,u), (14) 

where 1 is the identity matrix in flavor space. In contrast to the fully constrained NMSSM, the Higgs 
and singlet soft masses are obtained again by the tadpole equations, and also the NMSSM-specific 
trilinear soft couplings are treated as free parameters. Thus, necessary input parameters are 

mo, M 1 / 2 , tan/3, A 0 , A, k, A k , A x , p eff . (15) 


3. Spectrum generators for the NMSSM and their features 

3.1. Features 

The main features of the publicly available spectrum generators for the NMSSM are summarized 
in table [I] Here, two different kinds of spectrum generators are distinguished: dedicated tools which 

4 Note that, in eq. ]ip| , only one element of the matrix at a time can correspond to a pole Higgs mass, and 

only for an appropriate choice of p 2 . 

5 For further reference, we define also m? = m 2 , m 2 = in 2 - , in 2 = m 2 . 
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FlexibleSUSY 

NMSSMCALC 

NMSSMTools 

SOFTSUSY 

SPheno 

Code 

type 

using SARAH 

stand alone 

stand alone 

stand alone 

using SARAH 

language 

C++ 

Fortran 77 

Fortran 77 

C++ 

Fortran 90 



and 90 




Models 

No Z 3 

7 

X 

7 

7 

7 

GUT model 

/ 

X 

7 

7 

7 

CPV 

(/) 

/ 

X 

X 

(/) 

Thresholds 

scale(s) 

M z 

M t , Msusy 

M t , Msusy 

Mz 

Mz 

EW parame¬ 

full one-loop 

OS definitions 

full one-loop 

full one-loop 

full one-loop 

ters 






Yukawas 

full one-loop; 

one-loop 

one-loop 

full one-loop; 

full one-loop 


two-loop QCD 

(S)QCD; 

(S)QCD+Yukawa; two-loop 

two-loop QCD 



two-loop QCD 

two-loop QCD 

QCD; 






optionally 






two-loop 






SQCD 


strong gauge 

one-loop 

— 

— 

one-loop 

one-loop 


top+SUSY 



top+SUSY 

top+SUSY 



Higgs mass 

calculation 



scheme 

DR 

OS, DR 

DR 

DR 

DR 

one-loop 

full 

full 

full 

full 

full 

two-loop 

c+(<ab + at) 

QLsOt-t 

a s (ab + at) 

a s (atb + at) 

~ I - OLiOLj. 


+ MSSM 


+ MSSM 

+ MSSM 




SUSY 

masses 



one-loop 

/ 

X 

7 

7 

7 

momentum ef¬ 

/ 

X 

7 (a s ,a t ,ab 

7 

7 

fects 



only) 



Other observables 

decays 

X 

7 

via NMHDECAY 

via NMHDECAY 

7 

flavour, g — 2 

X 

X 

7 

X 

7 


Table 1: Public spectrum generators for the NMSSM and their features. The first two rows specify if the code is stand 
alone or gets necessary information from the Mathematica package SARAH, and in what language the code is written. The 
next three rows list the models that are supported: “GUT model” means that the code accepts as input a set of SUSY- 
breaking parameters at the GUT scale, and “CPV” shows if a Higgs sector with CP-even and -odd mixing is supported. 
Here, (S ) is used if this is only possible at the one-loop level so far. The rows under the sub-heading “Thresholds” give 
details of what SUSY and non-SUSY threshold corrections are included to calculate the running parameters entering the 
Higgs-mass calculation. The first row contains the scale (or scales) at which the threshold corrections are introduced, 
and the following three rows refer to different categories of parameters. Under “Higgs mass calculation” the main aspects 
of the calculation of the Higgs mass are given: the renormalization scheme (“scheme”), and what corrections at one- and 
two-loop are included, “full” at one-loop refers to a diagrammatic calculation including the entire p 2 dependence and all 
corrections up to those suppressed by first and second generation Yukawa couplings (which are also included in some 
codes). At two-loop the included corrections are shown. “+MSSM” refers to the approximation that the known MSSM 
results for the pure Yukawa interactions at two-loop are included. In the last column, ol% and aj stand for any of the 
couplings at, a 5 , a T , a\ and a K . Under “SUSY masses” it is shown if the code calculates also one-loop corrections to 
SUSY states (neutralinos, charginos, sfermions, gluino) and if the effects from the external momenta are included in these 
calculations. Finally, it is summarized in “Other observables” if the tool gives also a prediction for flavor observables 
and/or other precision observables like g — 2 and if it calculates the decays of at least the Higgs bosons. 
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contain all physically relevant information about the NMSSM (vertices, mass matrices, renormalization 
group equations (RGEs), loop corrections, etc.) in a hardcoded way, and multi-purpose generators 
which use the Mathematica package SARAH to derive this information for a given model. The advantage 
of the first approach, which is adopted by NMSSMTools, NMSSMCALC and SOFTSUSY, is that the tools 
can be used directly after compilation without further effort. In contrast, to use FlexibleSUSY and 
SPheno, one needs to first install SARAlf]] which is then used to generate the necessary source code to 
get an NMSSM spectrum generator. These additional steps allow these tools to be used for a much 
broader range of models, i.e. not only other singlet extensions beyond the NMSSM are supported, but 
also many other extensions of the MSSM. Also, the features of the NMSSM version (like boundary 
conditions at all relevant scales) can be adjusted in the input files for SARAH or FlexibleSUSY without 
the need to modify the final codcQ SARAH is available at hepforge 

http://sarah.hepforge.org/ 

More details on the different tools are given in the following, in the order in which the codes became 
publicly available. 

3.1.1. NMSSMTools 

The first public codes to be made available for the NMSSM were NMHDECAY [42] [43], NMSPEC [44j 
and NMGMSB, which are today collected in the package NMSSMTools. NMHDECAY was based on a suitably 
modified version of HDECAY [551 [54] [55] . This package is available at 

http://www.th.u-psud.fr/NMHDECAY/nmssmtools.html 

So far NMSSMTools is restricted to the real NMSSM with and without Z3 symmetry, but a version to 
support also CP violation is under construction |66| . NMSSMTools allows to calculate the Higgs masses 
for a parameter point with three different options: 

1. The most precise calculation makes use of the NMSSM corrections of ref. [36]. This provides 
a fully diagrammatic calculation of the Higgs masses at the one-loop level, including all con¬ 
tributions and the momentum dependence, plus the two-loop corrections of 0(a s (a b + at)) 
at zero external momentum. In addition, the two-loop corrections known from the MSSM at 
0((a t + a b ) 2 + a T (a T + a b )) [67116511661157] are included. 

2. The second option implements the same corrections as the first one, but neglects the momentum 
dependence of the one-loop self-energies. (This makes the code faster, which is useful for large 
scans.) 

3. The third option is the original implementation of the Higgs mass calculation in NMSSMTools, 
which - as described in the appendix C of ref. [5] - implements only a partial computation of the 
one-loop contributions, and restricts the two-loop contributions to the leading-log approximation. 

NMSSMTools also calculates the SUSY spectrum at the one-loop level, and the Higgs and sparticle 
decay branching fractions (the latter in NMSDECAY [70] using a suitable modification of SDECAY [71]). 
H-physics observables and the muon anomalous magnetic moment are computed following refs. [72173], 
and compared to experimental constraints. Bounds on the Higgs sector from LEP m and on the 
couplings of the SM-like Higgs boson (from ref. [75]) are implemented. A link to MicrOmegas is 
included which allows for an easy calculation of the dark matter relic density, and direct and indirect 
dark matter detection cross sections. 

It is possible to define a set of input parameters at the SUSY scale, or to define specific SUSY- 
breaking scenarios such as minimal supergravity using NMSPEC, or scenarios of gauge-mediated SUSY 
breaking (following ref. [76]) using NMGMSB. In all cases, the fine-tuning of the EWSB scale with respect 
to the input parameters can be estimated. 


^FlexibleSUSY also offers pre-generated code for selected models, with specific boundary conditions and EWSB 
outputs, allowing one to use these specific FlexibleSUSY spectrum generators without first installing Mathematica and 
SARAH. 

7 For easier comparison with other codes, we have created input files to allow for purely SLHA1 EU input for SPheno 
and FlexibleSUSY, while by default a block form for all matrices is used as in SLHA2 m- Interested readers can obtain 
the files from the authors. 
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3.1.2. SPheno and SARAH 

The main features of the SPheno version generated with SARAH for the NMSSM are a precise mass- 
spectrum calculation using the full two-loop RGEs of the NMSSM, with all flavor and CP effects, for 
the running between Mz, Msusy and Mqut- The RGEs are calculated using generic expressions 
given in Refs. [771 EH ESI [HOI ISH [HU ESI IM1- SPheno calculates also the full one-loop corrections to 
all masses, including the entire momentum dependence in the DR scheme. At two-loop all correc¬ 
tions in the gaugeless limit with vanishing external momenta are included [851|86| , by making use of 
generic expressions given in Refs. [571 [SHUSH]. This approach makes it possible to include all two- 
loop corrections without making use of any MSSM approximations as done by other codes, i.e. also 
the NMSSM-specific corrections of 0(a\(a\ + a K + a t )) are fully included. However, it is optionally 
also possible to make use of results from the literature for the two-loop corrections in the MSSM and 
NMSSM. These routines are usually much faster but have to be used carefully for models beyond the 
MSSM. Higgs sectors with CP violation are so far supported only at the one-loop level. 

Other features of the SPheno NMSSM version are the calculation of all important quark-flavor 
violating observables (B and K decays, AMs da /AAf^, e^') as well as lepton-flavor violating observ¬ 
ables using the FlavorKit functionality [50. Also, (g — 2)/ and electromagnetic dipole moments are 
predicted. Moreover, all two- and three-body decays of SUSY particles, and two-body decays of the 
Higgs scalars are calculated. The sparticle decays are purely tree-level, while for the Higgs bosons 
the NLO-QCD corrections to decays in two quarks, photons and gluons are included. Also decays in 
virtual vector bosons are taken into account. Finally, SPheno can make a prediction for the electroweak 
fine-tuning according to measures proposed in Refs. mm and writes all necessary input files to test 
points with HiggsBounds [S3] EH E5| an d HiggsSignals T)(i . 

The basic SPheno version that is used together with SARAH can be downloaded from 

http://spheno.hepforge.org/ 

A handy possibility to interface SPheno versions created by SARAH with Monte-Carlo tools like MadGraph 
HZ] EH] and WHIZARD [55-, lOOj as well as MicrOmegas is available via the BSM Toolbox [ 1011 . 

3.1.3. NMSSMCALC 

NMSSMCALC allows for the computation of the Higgs masses both in the CP-conserving and CP- 
violating NMSSM, both at one-loop and at two-loop level. In contrast to the other codes, NMSSMCALC 
makes use of mixed DR-OS renormalization conditions for the computation of the Higgs masses. In 
the (s)top sector the user can choose between OS and DR renormalization. In addition, if the charged 
Higgs boson mass Mjj± is given as an input instead of A\ , Mu± enters the Higgs mass calculation as 
an OS parameter. For comparison between the codes, however, the option of DR renormalization of 
the (s)top sector and the option with the input A\ are chosen. The Higgs mass calculation at one-loop 
level is performed including the full momentum dependence and all possible contributions [35] SE, ■ At 
the two-loop level the 0(a s at ) corrections are included [HU- They include the 0{a s at ) part relating 
the vacuum expectation value to physical observables, which is missing so far in the other spectrum 
generators. 

NMSSMCALC can be downloaded from 

http: //www. itp . kit. edu/^maggie/NMSSMCALC/ 

and comes together with an NMSSM extension of HDECAY |63l f64j [65] for the Higgs boson decays. 
The decays are calculated in the CP-conserving and CP-violating NMSSM. The decay widths include 
the dominant higher-order QCD corrections. Furthermore, the decays of the neutral Higgs bosons 
into a bottom pair include the higher order SUSY-QCD and the approximate SUSY-electroweak (EW) 
corrections up to one-loop order. The decays into a strange quark pair include the dominant resummed 
SUSY-QCD and the ones into r pairs the dominant resummed SUSY-EW corrections. Accordingly 
the charged Higgs boson decays into fermion pairs include the higher order SUSY corrections. All 
two-body decays into SUSY particles have been implemented. For the CP-conserving case the decays 
into stop and sbottom pairs come with the SUSY-QCD corrections. Finally, all relevant off-shell decays 
into two massive gauge boson final states, into gauge and Higgs boson final states, into Higgs pairs as 
well as into heavy quark pairs are included. For more details and a complete list of references see [4B,- 
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3.1.4■ SOFTSUSY 

SQFTSUSY is a widely used spectrum generator for the MSSM with and without i?-parity violation 
which has recently been extended to the NMSSM. The NMSSM spectrum generator is implemented to 
solve scenarios with and without the Z 3 symmetry. The implementation of EWSB conditions allows one 
to choose as output for the Z 3 -conserving case either the set { K,v s ,m .g} or the set {m 2 H . rn 2 H j . rri 2 }, 
and in the Z 3 -violating case there is an additional option of giving out {/ 1 , B^, % s }. For the high- 
scale boundary conditions there are a number of pre-defined options including the mSUGRA-inspired 
semi-constrained NMSSM and a general high-scale boundary condition which allows one to set all 
soft SUSY-breaking parameters independently. The users may also easily create their own boundary 
conditions. 

For the NMSSM, SOFTSUSY extends the full three-family one- and two-loop RGEs of the MSSM [771 
iina using general two-loop expressions from the literature [771 [75', with the NMSSM specialization 
given in Appendix D of the NMSSM manual m- These expressions were also cross-checked against 
two-loop NMSSM RGEs in the literature that use the 3rd-family approximation [9], and tested nu¬ 
merically against the RGEs from SARAH using an early development code from FlexibleSUSY. Finally, 
the RGEs for the vacuum expectation values are implemented using Refs. |831 [53]. The one-loop 
self-energies and tadpole corrections for the Higgs bosons were extended to the NMSSM using the 
expressions in Refs. EBI3Z|. The NMSSM extension also includes the two-loop corrections at order 
0(a s (a t + otb)) from Ref. |35| and in addition uses the MSSM results from Refs. [67,ESI 1S 3 IS?.. 

Higgs and sparticle decays may be obtained by interfacing SoftSUSY with NMHDECAY [42 43 and 
NMSDECAY uni, which is based upon Ref. m- These are both distributed as part of the NMSSMTools 
package, and SOFTSUSY provides a script to do the interface with this package automatically. 

The homepage of SOFTSUSY is 

http://softsusy.hepforge.org/ 


3.1.5. FlexibleSUSY 

FlexibleSUSY is a Mathematica and C++ package, which creates spectrum generators for user- 
defined (SUSY or non-SUSY) models. It makes use of SARAH to obtain model-dependent details, 
i.e. couplings, RGEs, self-energies and tadpole corrections. The boundary conditions on the model 
parameters at different scales, as well as spectrum-generator specific configuration details are specified 
in a FlexibleSUSY model file. FlexibleSUSY can therefore generate spectrum generators for the 
different MSSM singlet extensions provided by SARAH, for example the NMSSM with and without a Z 3 
symmetry. The generated spectrum generators consist of C++ classes, which can be easily adapted 
allowing one to use the routines in a custom-made program for solving non-standard problems. For 
example one may build a tower of effective field theories, which could be useful for NMSSM scenarios 
with very heavy superparticles. Currently, FlexibleSUSY fully supports only real parameters. Complex 
parameters have been developed very recently in version 1.1.0, but they are still undocumented and 
undergoing testing. The automatic implementation of alternative boundary-value-problem solvers and 
automatic generation of towers of effective field theories are currently under development. Additionally, 
an extension to calculate decays, FlexibleDecay, for both Higgs and superparticles is currently being 
developed. 

Due to the use of SARAH, FlexibleSUSY employs the full two-loop RGEs and one-loop corrections 
to the pole mass spectrum. At the two-loop level it makes use of Higgs mass corrections available in 
the literature. In the case of the NMSSM FlexibleSUSY uses the a s (a t + a^) corrections given in 
Ref. [36]. In addition, the MSSM two-loop corrections based on Refs. [571 EH [MUST] are used, though 
the user may choose to switch any of these two-loop corrections off using the SLHA file. 

The homepage of FlexibleSUSY is 


http://flexiblesusy.hepforge.org/ 





3.2. Used codes and options 

For the comparison in the following, these versions of the different codes were used: 

• NMSSMTools 4.7.0 

• SPheno 3.3.6 together with SARAH 4.5.7 

• NMSSMCALC 1.0.3 

• SOFTSUSY 3.6.0 

• FlexibleSUSY 1.1.0 together with SARAH 4.5.7 

We always make use of the most precise calculation available in each code, i.e. all available two-loop 
corrections turned on and in NMSSMTools we only compare Higgs masses calculated using option 1, 
described in section Em and do not consider any of the approximations from alternative options. 
We also concentrate on calculations done in the DR scheme. A discussion of the effects using an OS 
scheme is beyond the scope of this comparison and will be given elsewhere including also results from 
the upcoming NMSSM version of FeynHiggs [10311104 , 105111061 il(T7 i|. 

3.3. Predictions for the CP-even and CP-odd Higgs masses 

We discuss here the Higgs mass predictions of the various codes for six different test points (TP) 
with the following features: 

• MSSM-like point (TP1) 

• MSSM-like point with large stop splitting (TP2) 

• Point with light singlet and A close to the perturbativity limit (TP3) 

• Point with heavy singlet and A close to the perturbativity limit (TP4) 

• Point with slightly lighter singlet. Additional matter is needed for perturbativity; inspired by 
NMP9 HU] (TP5) 

• Point with huge A (TP6) 

In table [2] we list the input parameters that are different for the six benchmark points. As is usual 
in NMSSM studies, for the last four points we choose smallish values of tan/3, to avoid suppressing 
the NMSSM-specific contributions to the SM-like Higgs mass. For very large values of tan/3 additional 
differences between the codes are expected, since not all of them include the two-loop corrections 
involving the bottom and r Yukawa coupling. 

All values in this table apart from the one for tan /3 are defined at the scale Q. The shown value for 
tan/3 refers to the scale The remaining soft SUSY-breaking parameters, common to all points, 

are 


Mj = 1500 GeV , A T = 0 . (16) 

In addition, if possible via the input file, the codes were forced to calculate the masses directly at 
the input scale Q. This is possible for all tools except NMSSMTools. More details about the general 
handling of the scales in the different codes will be given in section |4.1.3| These points are chosen in a 
way that differences between the predicted Higgs masses are visible. However, they are not the result 
of an exhaustive scan to find points with the largest possible differences. Thus, we expect that these 
points are representative of a non-negligible fraction of the parameter space of the NMSSM. 


8 Since NMSSMCALC and NMSSMTools expect tan 0 to be defined at Q as well, the input values have to be adjusted 
accordingly. 


9 












Q 

tan j3 

A 

K 

Ha 

a k 

Me// 

Mi 

M 2 

m 3 

A-t 

Ab 

m i. 

m i„ 

TPl 

1500 

10 

0.1 

0.1 

-10 

-10 

900 

500 

1000 

3000 

3000 

0 

1500 

1500 

TP2 

1500 

10 

0.05 

0.1 

-200 

-200 

1500 

1000 

2000 

2500 

-2900 

0 

2500 

500 

TP3 

1000 

3 

0.67 

0.1 

650 

-10 

200 

200 

400 

2000 

1000 

1000 

1000 

1000 

TP4 

750 

2 

0.67 

0.2 

405 

0 

200 

120 

200 

1500 

1000 

1000 

750 

750 

TP5 

1500 

3 

0.67 

0.2 

570 

-25 

200 

135 

200 

1400 

0 

0 

1500 

1500 

TP6 

1500 

3 

1.6 

1.61 

375 

-1605 

614 

200 

400 

2000 

0 

0 

1500 

1500 


Table 2: Input parameters for the test points. All dimensionful parameters in units of GeV. 


TPl and TP2 have small couplings between the Higgs doublets and the singlet. Thus, a mass of 
about 125 GeV for the SM-like state can only be obtained via large (s)top corrections as in the MSSM. 
In contrast, TP3-TP6 come with a large A coupling and a small value for tan /3. This already increases 
the mass of the SM-like scalar at tree level. However, additional loop corrections beyond the MSSM 
are also expected to become more important. The difference between these four points is that TP3 
and TP4 have values for A and n that give a consistent model without any Landau pole up to the 
GUT scale of 10 16 GeV. For TP5 A becomes non-perturbative below the GUT scale unless additional 
matter is assumed to change the running. TP6 has a huge A which always leads to a cut-off. Scenarios 
with these large A couplings can be motivated by naturalness considerations [109 . 



TPl 

TP2 

TP3 

TP4 

TP5 

TP6 


hi 

FlexibleSUSY 

123.55 

122.84 

91.11 

127.62 

120.86 

126.46 

NMSSMCALC 

120.34 

118.57 

90.88 

126.37 

120.32 

123.45 

NMSSMTOOLS 

123.52 

121.83 

90.78 

127.30 

119.31 

126.63 

SOFTSUSY 

123.84 

123.08 

90.99 

127.52 

120.81 

126.67 

SPHENO 

124.84 

124.74 

89.54 

126.62 

119.11 

131.29 


h 2 

FlexibleSUSY 

1797.46 

5951.36 

126.58 

143.11 

125.08 

700.80 

NMSSMCALC 

1797.46 

5951.36 

124.86 

142.59 

123.14 

701.02 

NMSSMTOOLS 

1797.46 

5951.36 

127.28 

144.07 

126.95 

700.46 

SOFTSUSY 

1797.46 

5951.36 

126.59 

143.02 

125.12 

701.01 

SPHENO 

1798.01 

5951.35 

126.77 

144-04 

125.61 

689.30 


h'i 

FlexibleSUSY 

2758.96 

6372.08 

652.95 

467.80 

627.28 

1369.53 

NMSSMCALC 

2756.37 

6371.31 

652.48 

467.42 

627.00 

1368.68 

NMSSMTOOLS 

2758.51 

6345.72 

651.03 

466.38 

623.79 

1368.90 

SOFTSUSY 

2758.41 

6370.3 

652.78 

467.73 

627.14 

1369.19 

SPHENO 

2757.11 

6366.88 

651.21 

467.5 

624.02 

1363.02 


Table 3: Masses for the CP-even scalars (in GeV) for TP1-TP6 when using the spectrum generators “out-of-the-box”. 
The values correspond to the two-loop results obtained by the different tools. For NMSSMCALC the value using the DR 
scheme is shown and for NMSSMTools the value using the most precise calculation is given. The masses for the SM-like 
scalar are written in bold fonts, those for the singlet-like scalar in italics. 


The scalar masses calculated for these six points by the five spectrum generators as they come 
“out-of-the-box” are listed in table [3j and the absolute values of the entries of the rotation matrix Zu 
that correspond to the SM-like scalar are listed in tables |4][5j For TPl, TP2 and TP6 the singlet-like 
scalar is very heavy and decouples, for TP4 the singlet-like scalar is just slightly heavier than the 
SM-like scalar, while for TP3 and TP5 it is the lightest scalar. We find that, even for the MSSM-like 
points, the differences in the mass of the SM-like scalar can exceed 5 GeV. For the points with very 
large A the difference can be even larger. The masses for the other two scalars usually agree quite well 
between the codes. Only SPheno gives masses of the singlet-like states, for TP5 and TP6, which differ 
by 1-2% compared to all other tools. 
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\Zil\ 

TP1 

\Z i2 \ 

\z*\ 

\Zn\ 

TP2 

\z i2 \ 

\Z i3 i 

\z*\ 

TP3 

\Zi2\ 

\Za\ 

FlexibleSUSY 

0.1039 

0.9946 

0.0076 

0.1034 

0.9946 

0.0004 

0.2172 

0.1888 

0.9577 

NMSSMCALC 

0.1039 

0.9946 

0.0076 

0.1034 

0.9946 

0.0004 

0.2177 

0.1923 

0.9569 

NMSSMTOOLS 

0.1039 

0.9946 

0.0076 

0.1038 

0.9946 

0.0004 

0.2229 

0.2115 

0.9511 

SOFTSUSY 

0.1039 

0.9946 

0.0076 

0.1034 

0.9946 

0.0004 

0.2164 

0.1883 

0.958 

SPHENO 

0.104 

0.9945 

0.0074 

0.1035 

0.9946 

0.0004 

0.2269 

0.2265 

0.9472 

FlexibleSUSY 

0.0075 

0.0069 

0.9999 

0.0096 

0.0006 

1 . 

0.2814 

0.9274 

0.2466 

NMSSMCALC 

0.0075 

0.0068 

0.9999 

0.0095 

0.0006 

1 . 

0.2811 

0.9265 

0.2502 

NMSSMTOOLS 

0.0076 

0.0069 

0.9999 

0.0101 

0.0006 

0.9999 

0.2782 

0.9216 

0.2708 

SOFTSUSY 

0.0075 

0.0068 

0.9999 

0.0096 

0.0006 

1 . 

0.282 

0.9273 

0.246 

SPHENO 

0.0075 

0.0067 

0.9999 

0.0096 

0.0006 

1 . 

0.2725 

0.919 

0.2851 

FlexibleSUSY 

0.9946 

0.1039 

0.0068 

0.9946 

0.1034 

0.0096 

0.9347 

0.3231 

0.1483 

NMSSMCALC 

0.9946 

0.104 

0.0068 

0.9946 

0.1034 

0.0095 

0.9347 

0.3234 

0.1476 

NMSSMTOOLS 

0.9946 

0.104 

0.0068 

0.9945 

0.1038 

0.0102 

0.9338 

0.3255 

0.1484 

SOFTSUSY 

0.9946 

0.104 

0.0068 

0.9946 

0.1034 

0.0096 

0.9347 

0.3234 

0.1476 

SPHENO 

0.9945 

0.104 

0.0068 

0.9946 

0.1035 

0.0097 

0.935 

0.3228 

0.1468 


Table 4: Absolute values of the entries of the rotation matrix Zh corresponding to the SM-like scalar for TP1-TP3 
when using the spectrum generators “out-of-the-box”. The values correspond to the two-loop results obtained by the 
different tools. For NMSSMCALC the value using the DR scheme is shown and for NMSSMTools the value using the most 
precise calculation is given. 


4. Breaking down the differences 

f.l. Main differences 

The differences we have observed are larger than one usually finds for the MSSM when comparing 
the widely used tools that perform a DR calculation (571158]. We pinned down four sources which will 
be discussed in the following. The numerical impact of these differences is discussed in sec. |4.4| Further 
differences, such as e.g. whether first and second generation fermions are included into the one-loop 
computation (NMSSMCALC, SPheno, FlexibleSUSY) or not (NMSSMTools, SOFTSUSY), are numerically 
tiny. 

4-1.1. Renormalization scheme 

All codes but NMSSMCALC perform a pure DR renormalization. NMSSMCALC uses instead mixed DR- 
OS renormalization conditions at one-loop level. At two-loop level it offers an option to renormalize 
the (s)top sector not only OS but also DR. 

4-1-2. Calculation of the running parameters 

A very crucial point is how the DR parameters for the gauge and Yukawa couplings at the SUSY 
scale are calculated, since these couplings affect all loop calculations. 

Top Yukawa coupling. The most important parameter entering the loop calculations in the Higgs 
sector is the top Yukawa coupling. To calculate the value of the DR-renormalized Yukawa coupling at 
the SUSY scale, there are two different approaches implemented in the codes: 

1. SUSY thresholds at Mz- FlexibleSUSY, SOFTSUSY and SPheno adapt to the NMSSM the 
approach presented in Ref. m to calculate the running parameters of the MSSM Lagrangian at 
the SUSY scale. In this setup all one-loop SUSY and non-SUSY thresholds are calculated at the 
weak scale, identified with the zT-boson mass Mz, and also the known two-loop QCD corrections 
of (TTT1HT5| are included. The evolution of the parameters from Mz to the SUSY scale Msusy 
is done using two-loop SUSY RGEs. A detailed description of the calculation is given in Ref. f49j 
for FlexibleSUSY, Ref. [4Q for SOFTSUSY and in Ref. | I113 for SPheno. All three codes use in 
principle the same approach: 
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\Zil\ 

TP4 

\z i2 \ 

\Zisi 

\Zn\ 

TP5 

\Zi2\ 

\z i3 \ 

\z*\ 

TP6 

\z i2 \ 

\z i3 \ 

FlexibleSUSY 

0.4866 

0.7908 

0.3712 

0.2867 

0.4021 

0.8695 

0.2741 

0.9417 

0.195 

NMSSMCALC 

0.4869 

0.7849 

0.3832 

0.2967 

0.4433 

0.8459 

0.2746 

0.9416 

0.1947 

NMSSMTOOLS 

0.4855 

0.7495 

0.4501 

0.3072 

0.4711 

0.8269 

0.2743 

0.9416 

0.1953 

SOFTSUSY 

0.4867 

0.7902 

0.3725 

0.2863 

0.402 

0.8697 

0.2747 

0.9417 

0.1941 

SPHENO 

0.4835 

0.7646 

0.4262 

0.3033 

0.4703 

0.8288 

0.2712 

0.9407 

0.2039 

FlexibleSUSY 

0.0370 

0.4059 

0.9132 

0.2189 

0.8561 

0.4681 

0.2927 

0.1115 

0.9497 

NMSSMCALC 

0.0447 

0.4158 

0.9084 

0.2055 

0.8354 

0.5099 

0.2919 

0.1113 

0.9500 

NMSSMTOOLS 

0.0834 

0.4727 

0.8772 

0.1974 

0.8184 

0.5396 

0.2937 

0.1114 

0.9494 

SOFTSUSY 

0.0385 

0.4066 

0.9128 

0.2195 

0.8561 

0.4679 

0.2918 

0.1107 

0.9501 

SPHENO 

0.0682 

0.4525 

0.8892 

0.1958 

0.8204 

0.5372 

0.2945 

0.1206 

0.948 

FlexibleSUSY 

0.8728 

0.4581 

0.1683 

0.9327 

0.3245 

0.1575 

0.9161 

0.3174 

0.2451 

NMSSMCALC 

0.8723 

0.4594 

0.1674 

0.9326 

0.3251 

0.1568 

0.9162 

0.3177 

0.2443 

NMSSMTOOLS 

0.8703 

0.4634 

0.1670 

0.9309 

0.3291 

0.1584 

0.9157 

0.3177 

0.2460 

SOFTSUSY 

0.8727 

0.4587 

0.1675 

0.9327 

0.3249 

0.1568 

0.9162 

0.3176 

0.2444 

SPHENO 

0.8727 

0.4589 

0.1667 

0.9326 

0.3252 

0.1567 

0.9164 

0.3172 

0.2443 


Table 5: The same as Tab. [5] for TP4-TP6. 


the DR mass of the top quark is computed as 


DR _ vole 


+ rh t £(?h() 


(17) 


where E(p 2 ) is the (dimensionless) self-energy for momentum pj^] For other fermions, the relation 
becomes more complicated because of the tan (3 resummation which is performed by the codes. 
The DR Yukawa couplings are then computed at Mz by dividing the DR fermion masses by 
the appropriate VEVs. There is however a difference between the codes in the definition of the 
mass rht entering the correction in eq. (171: SOFTSUSY and FlexibleSUSY use to* = m ^ ole , while 
SPheno uses rht = mP R . The effect of the different definition of rht in the one-loop part of 
correction is partially compensated for by an appropriate shift in the two-loop QCD part of the 
correction. 


2. SUSY thresholds at Msusy- In this scheme, the pole top mass is converted to the MS 
scheme at the scale rrit, and the RGEs of the SM are used for the running to Msusy- At 
Msusy the conversion of the MS top mass of the SM to the corresponding DR Yukawa coupling 
of the NMSSM takes place including the SUSY thresholds. This is the approach followed in 
NMSSMTools and NMSSMCALC. However, there remain some differences in the concrete realization. 
NMSSMCALC uses the NLO-QCD contributions to the top-mass RGE for the running up to Msusy , 
and includes only the SQCD threshold corrections at Msusy- NMSSMTools resums also large 
logarithms ~ Y 2 log ( Msusy/ m t ) and ~ Y 2 log {Msusy /Ma) (the latter in case the scale Ma of 
the heavy MSSM-like scalars differs from the scale Msusy of the stop/sbottom) and it includes 
also O(at) threshold corrections from stop-Higgsino loops at Msusy- 


Strong interaction coupling. The strong interaction coupling at the scale Msusy , important for the 
two-loop corrections to the Higgs masses, is derived as follows: 


9 More precisely, SPheno calculates the chiral and scalar self-energy contributions Si, S fl , Eg to the 3x3 mass matrix, 
but this only leads to differences if flavor violation is included. 
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• FlexibleSUSY, SOFTSUSY, SPheno calculate a^ n (Mz) via 


,DR 


(Mz) = ^ 


(5),MS 


(Mz) 


A a.{M z y 


a tut \ a s ( 1 01 2 m t 1 

Ao,(M z ) =-I - - 21og — - s tog — - - X: 

\ i=l 


, ms, 

log Ife 


log w 


( 18 ) 

(19) 


where ai 5 ' ),MS denotes the strong coupling constant of the SM with Nf = 5 active flavors. Finally, 
a^(M z ) is run to Msusy using the two-loop RGEs of the NMSSM. 

• NMSSMCALC: o^ k (Msusy) is computed by running a^’ MS (Mz) up to the top pole mass m t 
using the QCD part of the two-loop RGEs of the SM with Nf = 5 active flavors. Afterwards, 
the so-obtained af s (mi) is run to Msusy using the corresponding RGEs with Nf = 6 active 
flavors. With this procedure the one-loop top threshold correction for a^ s is taken into account 
implicitly. Finally, a^ s (Msusy) is converted to DR by 


a» R (M SUS Y) = 


afHMsusv) 

_ pi™ 3 (Msusy) 

47r 


( 20 ) 


• NMSSMTools: a^ s ( Msusy ) is computed with aprocedure analogous to the one of NMSSMCALC, but 
it is not converted to the DR scheme. This difference amounts to a three-loop (i.e., higher-order) 
effect in the calculation of the Higgs masses. 

Electroweak interactions. The EW gauge sector of the NMSSM is determined by three fundamental 
parameters, i.e. the gauge coupling^ 10 ^ and g 2 and the combination of Higgs VEVs v = (v 2 + v^) 1 / 2 . 
For a consistent determination of these three parameters, a choice of three physical inputs is required. 
The procedures adopted by the five codes differ in the choices of physical inputs, as well as, in the case 
of NMSSMCALC, the choice of renormalization scheme for the fundamental parameters. 

1. FlexibleSUSY, SOFTSUSY and SPheno adapt again the approach of Ref. m to calculate g\ 
and g 2 . The three physical inputs are chosen as the Z mass, the Fermi constant Gf and the 
electromagnetic coupling of the SM at the scale M z in the 5-flavor scheme, ai 5 m’ MS (M z ). The 
threshold corrections to a ern from SUSY states, the charged Higgs and the top are calculated at 
the scale M z as 


<*(M Z ) = 


(5),MS 
Oi-em 


(M z ) 


1 - A a(M z y 


a /.an a (1 16, m ( 4 

A a (M z )=-(-^ y lo g—-- 


4 \ ' . ITlui 1 1 m di 


4 

3 


H lo S 


i—1 


m xt 

Mz 


1 

3 




*=l 


met 

Mz 



( 21 ) 


( 22 ) 


The Weinberg angle sin DR @w at the scale M z is obtained iteratively from the above-computed 
a^(M z ), plus Gp and M z , via the relation 



DR 


@w cos DR 0^j 


2 


k<£™(Mz) 

V2M 2 z G F (l-5 r )' 


( 23 ) 


10 For the U( 1) coupling gi we adopt the SM normalization, i.e. gi = g'. 
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where 


5 r p 


P = 


n ww (o) ^ e n l z (Ml) , jc i jc( 2) 

~T72 - 772 -+ d VB + <U b 




1-A/5’ 


Ml 


A p = 3?e 


Mzz(Mz) 

pMl 


n T 

11 ww 


(M^) 


Mir 


A p' 


,( 2 ) 


(24) 

(25) 


IIyy(p 2 ) (U = Z, W) are the DR-renormalized transverse parts of the self-energies of the 
vector bosons, computed at the renormalization scale Q = M z , and 5^ and Ap^ are two-loop 
corrections given in mm- The one-loop vertex and box corrections 5 vb implemented into 
FlexibleSUSY and SQFTSUSY are from refs. ma nm urn, generalized to give the complete 
one-loop result for the NMSSM, while the ones used in SPheno are automatically generated with 
SARAH for the NMSSM. We remark that, in this approach, the physical value of the W mass is 
an output of the calculation. At a generic scale Q, the relation between the pole mass My of a 
gauge boson V and its corresponding DR mass reads 

M™(Q) = yjMl + U^y(Ml) , (26) 

where the self-energy is also computed at the scale Q. 

In FlexibleSUSY sin DR Qw can be alternatively calculated from the DR W and Z masses, via 


sin Qw(Mz) = 


1 ( m™\m z ) 

\| \M^(M Z ) / 


(27) 


In this case the three physical inputs used to calculate g\ and g -2 are M z , M\y and ai 5 m’ MS (M z ). 
This approach is the original one implemented in FlexibleSUSY 1.0.0 and when we refer to 
“out-of-the-box” results for FlexibleSUSY it will be using this option. 

The so-obtained values for and sin DR 0^ are used to get g^ K (M z ) and g^ n (M z ), which 
are then evolved to the SUSY scale using two-loop RGEs of the NMSSM. 

2. NMSSMTools uses M z , M\y and Gf as physical inputs for the EW sector. It calculates the DR 
values of g\ and g 2 at the scale Msusy by making use of the relations between the fundamental 
EW parameters and the DR masses of the gauge bosons, the latter computed at the scale Q = 


Msusy as in eq. (261: 


V g? R (M SU SY) 2 + g^(MsusY) 2 =2 M% r (M S usy)/v dk (M susy ) , (28) 

g% R (M susy) = 2 M™(M S usy) / v dr (M S usy) , (29) 


with v db -(Msusy ) calculated via eq. (34 below. 

3. NMSSMCALC uses M z , My/ and the electromagnetic gauge coupling e as physical inputs for the 
EW sector. Moreover, it does not compute DR values for the EW gauge couplings, but adopts 
the OS conditions 

Mz 6 (30) 


9 1 = e 


My\r 


92 = 


V 1 ~ Mir,/Ml ’ 


where Myr, M z and e are the measured values. 


11 All expressions assume v ~ 246 GeV. 
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Electroweak VEV. The last SM parameter that enters the calculation of the Higgs masses is v, which 
gets translated into the two Higgs VEVs v,i and v u using the input value of tan /?. The parameter v at 
the scale Msusy is obtained as follows in the different codes: 

1. FlexibleSUSY and SOFTSUSY calculate the DR value at M z as 


v° R (M z ) 


2M^ r (M z ) 

y/gW(Mz) 2 + gW(Mz) 2 


(31) 


w DR (Msysy) is then obtained using two-loop RGEs with the ^-dependent terms calculated in 
the Feynman gauge [S3] lUj. 

2. SPheno calculates v DR directly at Msusy as 


v dr (M SU sy) 


2M rr (M SUS y) 

\Jg? R {M S usy ) 2 + g% R (MsusY ) 2 


(32) 


For completeness, we note that v° R (M z ), which is used in the threshold corrections to the gauge 
and Yukawa couplings, is calculated separately via the relation 


V DR (M Z ) = J M^ R (M Z ) 2 O_ ^ 2 e^) Sin20 ^ R ■ (33) 

In this way, v DR (M z ) and v° r (Msusy) are not strictly correlated via an RGE evolution, but 
differences of a higher order are present. 

3. NMSSMTools calculates v T:R (Msusy) via 

v dr (Msusy)~ 2 = V2G^ ^1 - W ]tf2 ^ - <>vb^ , (34) 

where <5vb denotes the one-loop correction to the muon decay amplitude. In the latter, NMSSMTools 
includes only the SM contribution. 

4. NMSSMCALC uses also for v an OS condition in terms of the three physical inputs: 

v = 2Mw yjl- Myy/Mz . (35) 


J^.1.3. Scales 

We have so far discussed differences in the threshold calculations, and how the codes perform those 
calculations at different scales. There are two other important scales involved in the calculation of the 
Higgs-boson masses: the scale at which the input parameters are taken by the codes, and the scale at 
which the pole masses are calculated. Thus, for a detailed comparison it is necessary that these scales 
be exactly the same in all codes. However, as mentioned in sec. |3.3[ it is not possible to unify these 
two scales among the codes without modifications in NMSSMTools. In the following, we describe how 
these two scales are treated in general by the different codes, and how - and to what extent - they 
can be fixed to particular numerical values using the provided options in the input. 

1. SPheno: If no scale is defined in the input, the input parameters are taken at Q = ^/ m i } m t 2 
and the masses are calculated at the same scale. If a scale is set via the entries EXTPAR[0] or 
M0DSEL[12], the input parameters are defined at this scale and the masses are calculated there 
as well. 

2. NMSSMCALC uses a similar approach to SPheno and always calculates the masses at the same scale 
where the input parameters are taken. However, if no explicit scale is defined, the default choice 

is Q = y/ m t L m iR . 
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3. SOFTSUSY by default takes the input scale to be the gauge-coupling unification scale. However, 
one may instead choose the input scale by either explicitly specifying it in EXTPAR[0] or by 
setting EXTPAR[0] to -1 in order to set the input scale to Q = \J m t^ m i 2 ■ I n contrast to SPheno, 
changing this entry does not affect the scale at which the pole masses are calculated. SOFTSUSY 
calculates the pole masses at Q = ^rnfjrjjT by default, like SPheno. To change this one may set 
an alternative scakp^ in the special SOFTSUSY block, using entry 4. 

4. FlexibleSUSY: Both the parameter input scale and the scale at which the pole masses are 
calculated may be chosen by the user in the model file used to generate the code. In the NMSSM 
model file used in this section, both scales are set via the entry EXTPAR[0]. In the model file 
used in sec. [5] to study GUT scenarios the parameter input scale is set to be the gauge-coupling 
unification scale, and the scale where the pole masses are calculated is set to a generalization of 
the SUSY scale which reduces to Q = ^/m^ rrif o in the absence of family mixing. 

5. NMSSMTools takes by default Q = yJ{2m\ L + m\ R + m 2 - )/2 for the input parameters and Q' = 

yjfni L mi R to calculate the Higgs, stop and sbottom pole masses. It is possible to change Q to 
a fixed value via EXTPAR[0] or MINPAR[0], but fixing Q' is only possible by modifying the code 
(which is done for the present study). The running between the two scales is performed with 
two-loop RGEs for the NMSSM. 

4-1.4- Included two-loop corrections 

While all codes use complete calculations of the one-loop corrections to the Higgs masses [55] [571 
15HJ1353, the partial two-loop corrections implemented in the codes cover different contributions. All 
the mentioned two-loop corrections are calculated for zero external momentum. 

1. NMSSMCALC includes only the two-loop corrections of 0(a s a t ) JH], but accounts also for the 
two-loop corrections to the electroweak VEV, which are omitted in the other codes. 

2. SOFTSUSY, FlexibleSUSY and NMSSMTools use the two-loop corrections of 0(a s (at -Pet;,)) which 
were calculated in the context of the NMSSM in Ref. [35 . In addition, these codes use the MSSM 
results for the two-loop corrections of 0((at+ab) 2 ) [55]. Moreover, SOFTSUSY and FlexibleSUSY 
include also routines to calculate the MSSM results for the tiny two-loop corrections involving 

a T |57] [53|55| [57] . 

3. SPheno performs a two-loop calculation in the gaugeless limit. Namely, any two-loop correction in 
the neutral Higgs sector that is independent of gi and (72 is included. However, SPheno includes 
also the available routines for the NMSSM corrections of 0(a s (at + ab)), as well as MSSM 
corrections involving only Yukawa couplings similar to SOFTSUSY and FlexibleSUSY. If wished 
by the user, these routines can be used instead of the full calculation. 


4-2. Adjusting the codes 

To check if the differences between the codes listed in the previous section can really explain the 
sizable discrepancies in the Higgs-mass predictions, we have modified the codes, choosing SOFTSUSY 
as reference. This choice resulted from mere convenience: FlexibleSUSY is very close to SOFTSUSY, 
and the differences between the two codes were already fully understood. This made it easy to adapt 
FlexibleSUSY and left just three codes which required more modifications. In particular, we have 
changed the following: 

• Calculation of the DR parameters: Changing the codes to mimic the calculation of the run¬ 
ning couplings of SOFTSUSY would be, of course, a big task beyond the scope of this comparisorp*] 
Therefore, the chosen approach was to read in directly the values of gi, Yi, v and tan/3 at the 
scale Q as calculated by SOFTSUSY and use them in the calculation of the Higgs masses. This 


12 For this scale to be used directly it must be larger than Mz- See documentation for full details. 

13 We stress that, in the context of this comparison, we double-checked the corresponding calculations performed by 
the codes and did not find any bugs. 
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original 

modified 


one-loop two-loop 

one-loop a s at a„(at + ctb) a s (at + a&) ‘full’ 

+ MSSM 

TPl 

FlexibleSUSY 

NMSSMCALC 

NMSSMTOOLS 

SOFTSUSY 

SPHENO 

119.48 123.55 

115.49 120.34 
123.52 

119.75 123.84 

120.69 124.84 

119.73 124.82 124.82 123.81 

119.75 124.85 

119.75 124.84 124.84 123.84 

119.75 124.84 124.84 123.84 

119.75 124.84 124.84 123.84 123.84 

TP2 

FlexibleSUSY 

NMSSMCALC 

NMSSMTOOLS 

SOFTSUSY 

SPHENO 

116.28 122.84 

111.80 118.57 

121.83 
116.49 123.08 

118.00 124.74 

116.46 123.79 123.79 123.06 

116.49 123.82 

116.49 123.82 123.82 123.08 

116.49 123.82 123.82 123.08 

116.49 123.81 123.81 123.08 123.05 

TP3 

FlexibleSUSY 

NMSSMCALC 

NMSSMTOOLS 

SOFTSUSY 

SPHENO 

124.16 126.58 

121.76 124.86 

127.28 
124.15 126.59 

124.79 126.77 

124.15 127.55 127.55 126.59 

124.15 127.56 

124.15 127.56 127.56 126.60 

124.15 127.56 127.56 126.59 

124.15 127.55 127.55 126.59 126.10 

TP4 

FlexibleSUSY 

NMSSMCALC 

NMSSMTOOLS 

SOFTSUSY 

SPHENO 

124.72 127.62 

122.85 126.37 

127.30 
124.59 127.52 

124.89 126.62 

124.59 128.23 128.23 127.51 

124.59 128.23 

124.59 128.24 128.24 127.52 

124.59 128.24 128.24 127.52 

124.59 128.23 128.23 127.52 126.33 

TP5 

FlexibleSUSY 

NMSSMCALC 

NMSSMTOOLS 

SOFTSUSY 

SPHENO 

122.53 125.08 

121.95 123.14 

126.95 

122.54 125.12 

122.90 125.61 

122.54 126.11 126.11 125.11 

122.54 126.12 

122.54 126.12 126.12 125.12 

122.54 126.12 126.12 125.12 

122.53 126.12 126.12 125.12 124.85 

TP6 

FlexibleSUSY 

NMSSMCALC 

NMSSMTOOLS 

SOFTSUSY 

SPHENO 

121.92 126.46 

118.55 123.45 

126.63 
122.12 126.67 

123.37 131.29 

122.09 127.7 127.7 126.64 

122.12 127.73 

122.12 127.73 127.73 126.67 

122.12 127.73 127.73 126.67 

122.09 127.71 127.71 126.65 129.91 


Table 6: The mass in GeV of the SM-like scalar after applying the modifications listed in sec. |4.2| for TP1-TP6 at 
the different loop levels and including different two-loop corrections. All two-loop corrections are for zero external 
momentum, while at one-loop the entire p 2 dependence is taken into account, “full” refers to a calculation which includes 
all corrections in the gaugeless limit where electroweak gauge couplings are neglected. Since the original version of 
NMSSMTools does not provide a flag to turn off the two-loop corrections, we don’t give any one-loop value here. 
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was done for NMSSMTools, NMSSMCALC and SPheno. In FlexibleSUSY no actual modifications to 
the code are made, but we will show results for both options for calculating the Weinberg angle. 
The option which matches the approach in SOFTSUSY will be treated as the modified version, 
while the approach implemented in FlexibleSUSY 1.0.0 will be treated as the original version. 

• Renormalization scale: NMSSMTools was forced to calculate the Higgs mass at the same scale 
as all other codes. 

• Renormalization scheme: The finite parts of the counterterms of the W boson mass SMw, 
the Z boson mass SMz, and the electric charge Se were put to zero in NMSSMCALC. Note that at 
two-loop level it is more convenient to introduce a counterterm for v, since in the gaugeless limit 
SM\ | and SM^ are zero, however, (5M§/M§ and SM^/Myy not. The remaining combinations 
can be expressed by Sv. This counterterm also needs to be set to zero to adjust to the other 
codes. 

• Loop corrections: Of course no additional loop corrections have been implemented in any of 
the codes in the context of this analysis. However, we added options to the different codes to 
turn off specific loop corrections, to get a clear picture of the masses at one-loop and at two-loop 
including (i) 0{a s a t ), (ii) 0{a s {a t +ab)) (iii) 0(a s (a t + «{,)) plus MSSM approximations for the 
corrections involving only Yukawa couplings, (iv) the full NMSSM calculation in the gaugeless 
limit. 

4-3. Masses after the adjustments 

The results for the SM-like Higgs mass after applying all changes to the codes are listed in Tab. [6] 
Here, we disentangled also the different loop corrections at the one- and two-loop level. The entries 
\Z([\ and \Z-$ \ of the corresponding rotation matrix are listed in Tab. [ 7 ] 

The big differences that were observed before are totally gone when comparing equivalent calcu¬ 
lations. Thus, for TP1-TP5 all codes show an impressive agreement at 0(a s ott). The corrections of 
0 (a s af) are completely negligible for all considered points (they could, however, give a sizable effect 
for very large tan 0). When including the purely Yukawa corrections at the two-loop level in the MSSM 
approximation, which is possible in all codes but NMSSMCALC, very good agreement is also found among 
the different tools. However, this approximation might not always be very good, and fails in particular 
for large A. This is manifest in the comparison with the full NMSSM result, which is available only 
for SPheno. 


4-4- Size of the different effects 

We discuss now in turn the impact of the different changes we applied. 

1. FlexibleSUSY: The entire difference between FlexibleSUSY and SOFTSUSY comes 
culation of the Weinberg angle sin DR Q\v at the scale Mz, as described in Section 


from 


4.1.2 


the cal- 


2. SPheno: In Tab.[8]the effects of adjusting successively the Yukawa couplings, the gauge couplings, 
and the electroweak VEV are shown. One can see that the biggest change comes from the Yukawa 
couplings. Of course, the top Yukawa coupling plays the main role. 

3. NMSSMCALC: Tab. [9] shows the effects of applying successively the different adjustments. Here 
we show the effect of changing the Yukawa couplings, changing the renormalization of the EW 
sector, and neglecting 0(a s at) contributions due to the conversion of the electroweak VEV to 
the DR parameter, which were omitted in the other spectrum calculators. The effect of the latter 
on the Higgs masses is proportional to A, and hence more visible for large A. 

4. NMSSMTools: Tab. [T0| shows the effects of applying successively the different adjustments: the 

effects of adjusting successively the scale Q' for the Higgs-mass calculation, the running Yukawa 
and gauge couplings and the running electroweak VEV. The changes in the scale have only an 
effect for TP2, because for the other points Q was chosen as which coincides with the 

scale Q' used by NMSSMTools. The only sizable effects are due to the top Yukawa coupling and 
the changes in the VEV. 
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original 

modified 


one-loop two-loop 

one-loop u s at a s (at + ctb) ‘full’ 

+ MSSM 

TPl 

FlexibleSUSY 

NMSSMCALC 

NMSSMTOOLS 

SOFTSUSY 

SPHENO 

0.1039,0.0076 0.1039,0.0076 

0.1039,0.0076 0.1039,0.0076 

0.1039,0.0076 
0.1039,0.0076 0.1039,0.0076 

0.1040,0.0076 0.1040,0.0074 

0.1039,0.0076 0.1039,0.0076 0.1039,0.0076 

0.1039,0.0076 0.1039,0.0076 

0.1039,0.0076 0.1039,0.0076 0.1039,0.0076 

0.1039,0.0076 0.1039,0.0076 0.1039,0.0076 

0.1039,0.0076 0.1039,0.0076 0.1039,0.0076 0.1039,0.0075 

TP2 

FlexibleSUSY 

NMSSMCALC 

NMSSMTOOLS 

SOFTSUSY 

SPHENO 

0.1034,0.0004 0.1034,0.0004 

0.1034,0.0004 0.1034,0.0004 

0.1038,0.0004 
0.1034,0.0004 0.1034,0.0004 

0.1035,0.0004 0.1035,0.0004 

0.1034,0.0004 0.1034,0.0004 0.1034,0.0004 

0.1034,0.0004 0.1034,0.0004 

0.1034,0.0004 0.1034,0.0004 0.1034,0.0004 

0.1034,0.0004 0.1034,0.0004 0.1034,0.0004 

0.1034,0.0004 0.1034,0.0004 0.1034,0.0004 0.1034,0.0004 

TP3 

FlexibleSUSY 

NMSSMCALC 

NMSSMTOOLS 

SOFTSUSY 

SPHENO 

0.2687,0.2975 0.2813,0.2470 

0.2664,0.3081 0.2811,0.2502 

0.2782,0.2708 
0.2693,0.2970 0.2820,0.2460 

0.2701,0.2951 0.2725,0.2851 

0.2693,0.2970 0.2835,0.2400 0.2819,0.2465 

0.2705,0.2935 0.2842,0.2385 

0.2705,0.2935 0.2842,0.2384 0.2828,0.2442 

0.2693,0.2970 0.2836,0.2399 0.2820,0.2460 

0.2693,0.2969 0.2836,0.2398 0.2820,0.2459 0.2719,0.2863 

TP4 

FlexibleSUSY 

NMSSMCALC 

NMSSMTOOLS 

SOFTSUSY 

SPHENO 

0.4852,0.3689 0.4866,0.3712 

0.4858,0.3625 0.4869,0.3832 

0.4855,0.4501 
0.4853,0.3700 0.4867,0.3725 

0.4852,0.3810 0.4835,0.4262 

0.4853,0.3698 0.4862,0.3938 0.4868,0.3722 

0.4845,0.4098 0.4836,0.4476 

0.4845,0.4098 0.4837,0.4472 0.4851,0.4235 

0.4853,0.3700 0.4862,0.3941 0.4867,0.3725 

0.4853,0.3698 0.4862,0.3939 0.4867,0.3723 0.4840,0.4133 

TP5 

FlexibleSUSY 

NMSSMCALC 

NMSSMTOOLS 

SOFTSUSY 

SPHENO 

0.0171,0.8787 0.2189,0.4681 

0.0154,0.9162 0.2055,0.5099 

0.1974,0.5396 
0.0208,0.8749 0.2195,0.4679 

0.0605,0.8192 0.1958,0.5372 

0.0205,0.8754 0.2445,0.3881 0.2191,0.4685 

0.1259,0.7020 0.2621,0.3276 

0.1259,0.7020 0.2621,0.3276 0.2485,0.3758 

0.0208,0.8749 0.2447,0.3876 0.2195,0.4679 

0.0206,0.8752 0.2449,0.387 0.2197,0.4672 0.1806,0.5766 

TP6 

FlexibleSUSY 

NMSSMCALC 

NMSSMTOOLS 

SOFTSUSY 

SPHENO 

0.2741,0.1946 0.2741,0.1950 

0.2742,0.1947 0.2746,0.1947 

0.2743,0.1953 
0.2746,0.1937 0.2747,0.1941 

0.2748,0.1938 0.2712,0.2039 

0.2746,0.1937 0.2747,0.1942 0.2747,0.1941 

0.2748,0.1935 0.2749,0.1939 

0.2748,0.1935 0.2749,0.1939 0.2749,0.1939 

0.2746,0.1937 0.2747,0.1941 0.2747,0.1941 

0.2746,0.1937 0.2747,0.1942 0.2747,0.1941 0.2711,0.2037 


Table 7: Absolute value of the Hd and singlet components (\Z^\, |Z^|) of the SM-like Higgs with index i after the 
modifications listed in sec. |4.2| for TP1-TP6. Because of the size of the table, we do not show explicitly the largest entry 
but it can easily be calculated as \Z^\ 2 = 1 — \Z^\ 2 — \Z ^| 2 . The remaining differences are caused by the external 
momentum p 2 used to calculate Z H : NMSSMCALC and NMSSMTools set p 2 = 0, while the other codes set p 2 = . The 

conventions for the given corrections are the same as in Tab.JSJ but we do not show the results for a s (a:b explicitly 

because they fully agree with the ones for a. s OLt for the number of digits used here. 


19 






























Point 

original 

Y 

9 

V 

modified 

TP1 

124.84 

123.65 

123.61 

123.84 

123.84 

TP2 

124.74 

123.18 

123.13 

123.05 

123.05 

TP3 

126.77 

126.06 

126.00 

126.10 

126.10 

TP4 

126.62 

126.21 

126.16 

126.33 

126.33 

TPS 

125.61 

124.89 

124.84 

124.85 

124.85 

TP6 

131.29 

130.06 

130.01 

129.91 

129.91 


Table 8: The Higgs prediction for the SM-like Higgs mass by SPheno after applying successively the different adjustments 
for the Yukawas ( Y ), gauge couplings (g). and the electroweak VEV (u). Here, “original” refers to the results when using 
the code without any modification, while for “modified” all adjustments are turned on. 


Point 

original 

Y 

6 i 

d 2 

modified 

TP1 

120.34 

124.41 

124.85 

124.85 

124.85 

TP2 

118.57 

123.31 

123.82 

123.82 

123.82 

TP3 

124.86 

127.55 

127.50 

127.56 

127.56 

TP4 

126.37 

128.32 

128.18 

128.23 

128.23 

TPS 

123.14 

126.21 

126.03 

126.12 

126.12 

TP6 

123.45 

127.26 

127.55 

127.73 

127.73 


Table 9: The prediction for the SM-like Higgs mass by NMSSMCALC after adjusting the DR Yukawa couplings (Y); 
performing a DR renormalization of the EW sector by changing the values of gi, <72 and v and removing the finite parts 
of the one-loop counter-terms 5M\y, 8Mz and <5e (<5i); removing also the finite counterterm for v at the two-loop level 
(<$ 2 ). The same conventions as for Tab.[8]are used. 


One can see from this comparison that the main source of discrepancies between the predictions 
of the five codes is the DR-renormalized top Yukawa coupling used in the calculation of the Higgs 
masses: differences in the determination of Y t DR (Q) have a larger impact than additional two-loop 
corrections to the Higgs mass matrix. In particular, the fact that NMSSMCALC includes only (S)QCD 
corrections in the calculation of Y t DR (Q) results in a prediction for the SM-like Higgs mass a few 
GeV below the ones of the other codes. The results in tables [6] and [9] suggest that the inclusion of 
additional one-loop corrections in the calculation of that coupling by NMSSMCALC would account for the 
bulk of the corresponding two-loop corrections to the Higgs masses, thus reducing the discrepancies 
with the other codes. However, the precise calculation of the Higgs decays in NMSSMCALC includes only 
(S)QCD corrections. Therefore, with the aim of preserving the overall consistency of the calculation, 
the authors of NMSSMCALC include only (S)QCD corrections in Y t DR (Q) as well. 


Point 

original 

<q 

ii 

<e> 

Y 

9 

V 

modified 

TP1 

123.52 

123.52 

123.96 

123.99 

123.84 

123.84 

TP2 

121.83 

121.44 

123.46 

123.52 

123.08 

123.08 

TP3 

127.28 

127.28 

127.43 

127.43 

126.60 

126.60 

TP4 

127.30 

127.30 

127.13 

127.07 

127.52 

127.52 

TPS 

126.95 

126.95 

127.34 

127.45 

125.12 

125.12 

TP6 

126.63 

126.63 

127.56 

127.66 

126.67 

126.67 


Table 10: Changes in the prediction by NMSSMTools for the SM-like Higgs mass after forcing the SUSY scale and the 
scale for the mass calculation to be identical ( Q' = Q ), changing the Yukawa couplings (Y), the gauge couplings ( g ), 
and the electroweak VEV ( v ). The same conventions as for Tab. [8] are used. 
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5. Expected effects for GUT scenarios 


We have focused so far on the Higgs-mass calculation in which all SUSY parameters are fixed by 
the input at the SUSY scale. However, most codes offer also the possibility to study a GUT scenario 
where the soft SUSY-breaking parameters are fixed at the GUT scale. The differences among the codes 
discussed in the previous sections are going to play an important role here as well. Indeed, because of 
the additional RGE running, the discrepancies between the results of the different codes are usually 
of the same size as, and often even bigger than those found in the case of SUSY-scale inputs with 
comparable parameter choices. In particular, the differences in the determination of the DR couplings 
can play a crucial role, because the different values for the couplings at the SUSY scale affect the entire 
running between M$usy and Mqut ■ This is especially important for points in which the input value 
of A at MsuSy is close to the perturbativity limit. Small changes in the threshold corrections to the 
weak-scale couplings can change the value of A at the GUT scale considerably. This will then change 
the running of all soft terms from the GUT scale to the SUSY scale. While the modifications in the 
sfermion parameters are not necessarily large, A\ and A K can be strongly affected. This will then have 
an impact on the singlet-like scalar. 

To demonstrate this effect, we consider the two points given in Tab. [Tl] Both points provide 
an SM-like Higgs mass in the preferred range together with a light singlet, but the first (GTP1) is 
characterized by small values of A and k, while the second (GTP2) is characterized by larger values 
of those couplings, not far below the perturbativity limit. To show the large sensitivity of GTP2 
compared to GTP1 to the parameters at the SUSY scale, we vary the top mass within the uncertainty: 
rn t = 173.2 ± 1.5 GeV. The masses of the three CP-even Higgs bosons computed by the four codes that 
allow for GUT-scale boundary conditions, i.e. FlexibleSUSY, NMSSMTools, SOFTSUSY and SPheno, are 
shown in Tab. [12] We have used for this comparison the un-tuned versions of all codes. 



m 0 

M\/2 

Ao 

tan/3 

A 

hi 

A\ 

a k 

Meff 

GTP1 

1000 

1000 

-3000 

10 

0.05 

0.05 

-150 

-300 

100 

GTP2 

1000 

1000 

-3000 

10 

0.48 

0.46 

-150 

-300 

100 


Table 11: Input values for GTP1 and GTP2. All dimensionful parameters in units of GeV. The parameters in the left 
section of the table are defined at the GUT scale, whereas those in the right section are defined at the scale MguSYi 
with the exception of tan [3 which is defined at Mz • 



rn t = 171.7 GeV m t = 173.2 GeV m t = 174.7 GeV 

FlexibleSUSY 

NMSSMTools 

SOFTSUSY 

SPheno 

100.1, 122.2, 951.9 100.0, 123.3, 967.0 100.0, 124.5, 982.0 

100.0, 119.3, 982.1 100.0, 120.1, 974.4 99.9, 120.8, 1014.0 

100.1, 122.5, 958.6 100.0, 123.6, 974.4 100.0, 124.8, 990.2 

100.0, 123.7, 976.6 99.9, 124.9, 992.3 99.9, 126.0, 1008.0 


m t = 171.7 GeV m t = 173.2 GeV m t = 174.7 GeV 

FlexibleSUSY 

NMSSMTools 

SOFTSUSY 

SPheno 

102.8, 123.3, 961.9 97.9, 124.6, 974.1 92.5, 126.1, 986.1 

96.6, 118.8, 987.7 91.4, 119.1, 1000.8 no output 

102.4, 123.7, 968.2 97.3, 125.2, 981.1 91.4, 126.7, 994.0 
95.1, 125.3, 981.9 89.3, 126.8, 994.6 82.5, 128.3, 1007.5 


Table 12: The masses of the three CP-even Higgs bosons, mj 1/1 2 3 , for GTP1 (top) and GTP2 (bottom) for three different 
values of the top quark mass, m*. 


We see that, for GTP1, the SM-like Higgs mass varies by about 2 GeV for all codes when mt 
is varied within one standard deviation, while the light-singlet mass is fairly independent of the top 
mass, and the predictions of all codes for the light-singlet mass agree within 0.1 GeV. For GTP2, the 
variation with m t in the SM-like Higgs mass is about 3 GeV, but the singlet is very sensitive to the 
value of the top mass: its mass can change by about 10 GeV, and the different codes predict masses 
that differ from each other by several GeV. In addition, for mt = 174.7 GeV NMSSMTools runs into 
the Landau pole and is not able to calculate the spectrum. The situation can be even more dramatic 
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when taking for instance the benchmark points proposed in Ref. m ■■ for each of the three points 
BP1-BP3, at least two codes do not produce an output. When they are forced to produce an output 
(for instance by ignoring tachyonic states which appear at tree-level), the differences between the codes 
can reach 20-30 GeV for the light singlet. 

Points with very large A at the SUSY scale are attractive since they seem to soften the SUSY fine- 
tuning problem (TTJ Q2J ED HU [U3 ED d3 ED ED HO] ■ However, if A is too close to the perturbativity 
limit, the output produced by the spectrum generators has to be taken with some care: these points 
can be very sensitive to small variations of the input parameters. Changing the SM input parameters 
within their respective uncertainties - or changing the details of the determination of the running 
couplings - can lead to a completely different spectrum, or even result in the spectrum generator 
classifying the parameter point as unphysical. 


6. Conclusions 

We have discussed the differences in the predictions for the scalar Higgs masses and mixings by 
public spectrum generators for the NMSSM. It has been shown that sizeable discrepancies can show 
up for various parameter points, but the origins of those discrepancies are fully understood. The main 
sources of discrepancies are in the approach used to calculate the running parameters that enter the 
Higgs-mass calculation. The most important parameter here is the top Yukawa coupling. Further 
discrepancies arise from the different accuracy of the two-loop corrections included in the codes. 

This has important consequences for phenomenological studies: For MSSM-like points, FlexibleSUSY, 
NMSSMTools, SOFTSUSY and SPheno usually agree quite well with each other, while NMSSMCALC returns 
in general a smaller value for the SM-like Higgs boson mass. The main reasons for this discrepancy are 
the value of Y t DR at Msusy , which is always smaller in NMSSMCALC, and the inclusion of the MSSM 
results for the 0((a t + cp,) 2 ) corrections in FlexibleSUSY, NMSSMTools and SOFTSUSY but not in 
NMSSMCALC. For MSSM-like points, the MSSM results provide a good approximation to the calculation 
implemented in SPheno which includes the NMSSM specific effects. On the other hand, this approxi¬ 
mation usually fails for large values of A and/or if light singlet states are present. In these cases, the 
results of FlexibleSUSY, NMSSMTools and SOFTSUSY have to be taken with care because it might be 
that the MSSM approximation has even the wrong sign compared to the full NMSSM calculation [40.. 
For these cases SPheno, with the additional NMSSM specific contributions, is expected to provide the 
most reliable prediction, with the caveat that if states in the loop become very light, the calculation 
can be plagued by intrinsic problems of the effective-potential calculation [1191 112011121 . NMSSMCALC, 
on the other hand, is the only code that includes two-loop 0 {atOt s ) corrections to v, which become 
non-negligible for large values of A. 

It has been demonstrated with one example, that the differences can be much more pronounced 
for GUT scenarios with values of A close to the perturbativity limit. The differences in the threshold 
calculations among the codes can cause huge discrepancies in particular for light singlets and the 
predicted masses can differ by 10 GeV and more. 

In general, comparing the predictions of different codes may provide a ballpark indication for the 
theoretical uncertainty of the Higgs mass calculation. However, even if these codes agree very well, one 
should not necessarily assume that the uncertainty is small: it is known from the MSSM that spectrum 
generators performing a DR calculation (SOFTSUSY, SPheno, Suspect [ 122] I can agree quite well, while 
sizeable differences to the OS calculation of FeynHiggs exists. The differences are assumed to come 
from the missing electroweak corrections and momentum dependence at two-loop level as well as from 
the dominant three-loop corrections. These are the effects that underlie the often-quoted estimate 
of a 3 GeV uncertainty for the SM-like Higgs mass in the MSSM [571 |1~(J5] . To get an estimate of 
the remaining theoretical uncertainty in the NMSSM, a comparison between different renormalization 
schemes will therefore be necessary. This is deferred to future work, when further publicly available 
codes will have implemented OS renormalization. However, even the validity of such an estimate is 
under discussion already in the MSSM, in particular for a heavier SUSY scale of 2 TeV and larger. 
Calculations based on the effective field theory approach, which become valid in this mass range, 
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usually predict a SM-like Higgs mass that is lower by a few GeV compared to the prediction from 
dedicated tools for the MSSM |T23 , 15311125 ]. 

With this work an important first step has been taken to understand the large differences in 
the Higgs mass computations of the various public NMSSM codes. For the investigated codes these 
differences are now completely understood and it has been shown that discrepancies of more than 
3 GeV are not unusual taking into account the different corrections implemented. To further pin down 
the residual theoretical uncertainty the next natural step will be to extend this investigation to the 
comparison of different renormalization schemes. 

Acknowledgements 

We thank Peter Drechsel, Mark D. Goodsell, Werner Porod, Heidi Rzehak, Michael Spira, Kathrin 
Walz and Georg Weiglein for many useful discussions. 

This work was supported in part by the Research Executive Agency (REA) of the European Com¬ 
mission under the Initial Training Networks “INVISIBLES” (PITN-GA-2011-289442) and “Higgs-Tools” 
(PITN-GA-2012-316704), and by the European Research Council (ERC) under the Advanced Grant 
“Higgs@LHC” (ERC-2012-ADG_20120216-321133). The work of P. S. at LPTHE is supported in 
part by French state funds managed by the ANR (ANR-ll-IDEX-0004-02) in the context of the ILP 
LABEX (ANR-10-LABX-63). The work of P. A. was supported by the ARC Centre of Excellence for 
Particle Physics at the Terascale. 


Appendix A. Used options in input files 

Appendix A.l. FlexibleSUSY 


Block 

MODSEL 

# 

Select model 

6 

0 

# 

flavour violation 

12 

SCALE 

# 


Block 

FlexibleSUSY 



0 

1.000000000e-04 

# 

precision goal 

1 

0 

# 

max. iterations (0 = automatic) 

2 

0 

# 

algorithm (0 = two_scale, 1 = lattice) 

3 

0 

# 

calculate SM pole masses 

4 

2 

# 

pole mass loop order 

5 

2 

# 

EWSB loop order 

6 

2 

# 

beta - functions loop order 

7 

1 

# 

threshold corrections loop order 

8 

1 

# 

Higgs 2-loop corrections 0 (alpha_t ^ 


alpha_ s ) 



9 

1 

# 

Higgs 2-loop corrections 0(alpha_b ^ 


alpha_ s ) 



10 

1 

# 

Higgs 2-loop corrections 0(alpha_t~2 


+ alpha_t alpha_b 

+ alpha_b ~2) 

11 

1 

# 

Higgs 2-loop corrections 0(alpha_tau~2) 


Appendix A.2. NMSSMCALC 

Block MODSEL 

31 # NMSSM 

5 0 # 0: CP - conserving; 2: general 

6 2 # loop level 1: one 2: two 

71 # 1: DRbar scheme for top/stop 

top/stop-sector 


CP -vio l ation 

sector; 2: OS scheme for 
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Appendix A. 3. NMSSMTools 


BLOCK MODSEL 

31 # NMSSM particle content 

10 # IMOD (0 = g eneral NMSSM, l=mSUGRA, 

2=GMSE) 

10 0 # ISCAN (0=no scan, l=grid scan, 

2=random scan, 3=MCMC) 

90 # IOMGFLAG1=0: no (default), =1: 

relic density only, 

# =2: dir. det . rate, =3: indir. 
det . rate, =/: both, 

# OMGFLAG >0 : 0 . 1 07<0MG <0.131, <0: ^ 

OMG<0.131 

8 2 # Precision for Higgs masses 

(default 0: as before, 

# 1: full 1 loop + full 2 loop from 


top/bot Yukawas 

courtesy of P. Slavich) 

13 0 
(default 0) 

11 0 

0=no, default=1) 

14 0 

H->VV, VV*, V*V* 

15 0 


BLOCK MINPAR 
0 


2: as 1 + pole masses - 1&2 by 


# 1: Sparticle decays via NMSDECAY 


# Constraints on (g-2)_muon (l=yes, 

# 0: H-> VV, VV * (default); 1: ^ 

# Precision for micromegas (defalt=0): 
+0/1: fast computation on/off 

+0/2: Beps=ld-3, Id-6 
+ 0//: virtual lVs off/on 


SCALE # MSUSY (If =/= SQRT (2*MQ1+MU 1+MD1)/2) 


Appendix A./. SOFTSUSY 


# Example NMSSM input in SLHA2 format 


Block MODSEL 
6 0 

1 0 

3 1 

12 SCALE 

Block SOFTSUSY 

1 1.000000000e-04 

2 2.000000000e+00 

5 1.000000000E+00 

3 1.000000000E+00 
4 SCALE 

7 4.0 

15 1.000000000E+00 

(defauIt: 0) 


# Select model 

# flavour violation 

# mSUGRA 

# NMSSM 

# 

# SOFTSUSY specific inputs 

# tolerance 

# up-quark mixing (=1) or down (=2) 

# 2-loop running 

# printout 

# 

# NMSSMTools compatible output 







15 


16 

17 

18 

19 

20 


21 


16 0.000000000E+00 


17 1.000000000E+00 

(defauIt: 0) 

18 1.000000000E+00 


# Call micrOmegas (default: 0 = no, 

# 1 = relic density only, 

# 2 = direct detection + relic density, 

# 3 = indirect detection + relic density, 

# 4 = all) 

# Sparticle decays via NMSDECAY 

# use soft Higgs masses as EWSB output 


Appendix A. 5. SPheno 


1 

Block MODSEL 


# 



2 

1 0 



# 


1/0: High/low scale input 

3 

2 1 



# 

Boundary Condition 

4 

6 0 



# 

Generation Mixing 

5 

12 

SCALE 

# 




6 

Block SPhenoInput 


# 

SPheno specific input 

7 

1 

-1 



# 

error level 

8 

2 

0 



# 

SPA conventions 

9 

7 

0 



# 

Skip 2-loop Higgs corrections 

10 

8 

3 



# 

Method used for two-loop calculation 

11 

9 

1 



# 

Gaugeless limit used at two-loop 

12 

10 

0 



# 

safe - mode used at two-loop 

13 

11 

0 



# 

calculate branching ratios 

14 

13 

1 



# 

include 3-Body decays 

15 

12 

1.000E-04 



# 

write only branching ratios larger than this “4 



value 





16 

15 

1.000E-30 



# 

write only decay if width larger than this value 

17 

31 

-1 



# 

fixed GUT scale (-1: dynamical GUT scale) 

18 

32 

0 



# 

Strict unification 

19 

34 

1.000E-04 



# 

Precision of mass calculation 

20 

35 

80 



# 

Maximal number of iterations 

21 

36 

40 



# 

Minimal number of iterations before 



discardi 

ng p 

oint 


22 

37 

1 



# 

Set Yukawa scheme 

23 

38 

2 



# 

1- or 2-Loop RGEs 

24 

50 

1 



# 

Majorana phases: use only positive masses 

25 

51 

0 



# 

Write Output in CKM basis 

26 

52 

1 



# 

Write spectrum in case of tachyonic states 

27 

55 

1 



# 

Calculate one loop masses 

28 

57 

0 



# 

Calculate low energy constraints 

29 

65 

1 



# 

Solution tadpole equation 


References 

[1] S. Chatrchyan et al. . “Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC,” 
Phys.Lett. , vol. B716, pp. 30-61, 2012, 1207.7235. 

[2] G. Aad et al., “Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS 
detector at the LHC,” Phys.Lett., vol. B716, pp. 1-29, 2012, 1207.7214. 

[3] J. Camargo-Molina, B. O’Leary, W. Porod, and F. Staub, “Stability of the CMSSM against sfermion VEVs,” 
JHEP , vol. 1312, p. 103, 2013, 1309.7212. 

[4] N. Blinov and D. E. Morrissey, “Vacuum Stability and the MSSM Higgs Mass,” JHEP, vol. 1403, p. 106, 2014, 
1310.4174. 


25 




[5] D. Chowdhury, R. M. Godbole, K. A. Mohan, and S. K. Vempati, “Charge and Color Breaking Constraints in 
MSSM after the Higgs Discovery at LHC,” JHEP, vol. 1402, p. 110, 2014, 1310.1932. 

[6] J. Camargo-Molina, B. Garbrecht, B. O’Leary, W. Porod, and F. Staub, “Constraining the Natural MSSM through 
tunneling to color-breaking vacua at zero and non-zero temperature,” Phys.Lett., vol. B737, pp. 156-161, 2014, 
1405.7376. 

[7] U. Chattopadhyay and A. Dey, “Exploring MSSM for Charge and Color Breaking and Other Constraints in the 
Context of Higgs@125 GeV,” JHEP , vol. 1411, p. 161, 2014, 1409.0611. 

[8] M. Maniatis, “The Next-to-Minimal Supersymmetric extension of the Standard Model reviewed,” Int. J. Mod. 
Phys. , vol. A25, pp. 3505-3602, 2010, 0906.0777. 

[9] U. Ellwanger, C. Hugonie, and A. M. Teixeira, “The Next-to-Minimal Supersymmetric Standard Model,” 
Phys.Rept. , vol. 496, pp. 1-77, 2010, 0910.1785. 

[10] U. Ellwanger and C. Hugonie, “The Upper bound on the lightest Higgs mass in the NMSSM revisited,” 
Mod.Phys.Lett., vol. A22, pp. 1581-1590, 2007, hep-ph/0612133 

[11] M. Bastero-Gil, C. Hugonie, S. King, D. Roy, and S. Vempati, “Does LEP prefer the NMSSM?,” Phys.Lett. , 
vol. B489, pp. 359-366, 2000, hep-ph/0006198 

[12] R. Dermisek and J. F. Gunion, “Consistency of LEP event excesses with an h -> aa decay scenario and low-fine- 
tuning NMSSM models,” Phys.Rev. , vol. D73, p. 111701, 2006, hep-ph/0510322 

[13] U. Ellwanger, G. Espitalier-Noel, and C. Hugonie, “Naturalness and Fine Tuning in the NMSSM: Implications of 
Early LHC Results,” JHEP , vol. 1109, p. 105, 2011, 1107.2472. 

[14] G. G. Ross and K. Schmidt-Hoberg, “The Fine-Tuning of the Generalised NMSSM,” Nucl.Phys., vol. B862, pp. 710- 
719, 2012, 1108.1284. 

[15] G. G. Ross, K. Schmidt-Hoberg, and F. Staub, “The Generalised NMSSM at One Loop: Fine Tuning and Phe¬ 
nomenology,” JHEP , vol. 1208, p. 074, 2012, 1205.1509. 

[16] T. Gherghetta, B. von Harling, A. D. Medina, and M. A. Schmidt, “The Scale-Invariant NMSSM and the 126 GeV 
Higgs Boson,” JHEP , vol. 1302, p. 032, 2013, 1212.5243. 

[17] M. Perelstein and B. Shakya, “XENON100 implications for naturalness in the MSSM, NMSSM, and A- 
supersymmetry model,” Phys.Rev., vol. D88, no. 7, p. 075003, 2013, 1208.0833. 

[18] D. Kim, P. Athron, C. Balazs, B. Farmer, and E. Hutchison, “Bayesian naturalness of the CMSSM and CNMSSM,” 
Phys.Rev. , vol. D90, no. 5, p. 055008, 2014, 1312.4150. 

[19] A. Kaminska, G. G. Ross, K. Schmidt-Hoberg, and F. Staub, “A precision study of the fine tuning in the DiracN- 
MSSM,” JHEP , vol. 1406, p. 153, 2014, 1401.1816. 

[20] M. Y. Binjonaid and S. F. King, “Naturalness of scale-invariant NMSSMs with and without extra matter,” 
Phys.Rev. , vol. D90, no. 5, p. 055020, 2014, 1403.2088. 

[21] H. Dreiner, F. Staub, and A. Vicente, “General NMSSM signatures at the LHC,” Phys.Rev., vol. D87, no. 3, 
p. 035009, 2013, 1211.6987. 

[22] S. King, M. Muhlleitner, R. Nevzorov, and K. Walz, “Natural NMSSM Higgs Bosons,” Nucl.Phys., vol. B870, 
pp. 323-352, 2013, 1211.5074. 

[23] S. King, M. Muhlleitner, R. Nevzorov, and K. Walz, “Discovery Prospects for NMSSM Higgs Bosons at the 
High-Energy Large Hadron Collider,” Phys.Rev. , vol. D90, no. 9, p. 095014, 2014, 1408.1120. 

[24] U. Ellwanger and A. M. Teixeira, “NMSSM with a singlino LSP: possible challenges for searches for supersymmetry 
at the LHC,” JHEP , vol. 1410, p. 113, 2014, 1406.7221. 

[25] J. Espinosa and M. Quiros, “Two loop radiative corrections to the mass of the lightest Higgs boson in supersym¬ 
metric standard models,” Phys.Lett., vol. B266, pp. 389-396, 1991. 

[26] J. Espinosa and M. Quiros, “Upper bounds on the lightest Higgs boson mass in general supersymmetric Standard 
Models,” Phys.Lett., vol. B302, pp. 51-58, 1993, hep-ph/9212305 

[27] T. Elliott, S. King, and P. White, “Supersymmetric Higgs bosons at the limit,” Phys.Lett. , vol. B305, pp. 71-77, 
1993, hep-ph/9302202 

[28] U. Ellwanger, “Radiative corrections to the neutral Higgs spectrum in supersymmetry with a gauge singlet,” 
Phys.Lett. , vol. B303, pp. 271-276, 1993, hep-ph/9302224 

[29] P. Pandita, “Radiative corrections to the scalar Higgs masses in a nonminimal supersymmetric Standard Model,” 
Z.Phys. , vol. C59, pp. 575-584, 1993. 

[30] T. Elliott, S. King, and P. White, “Squark contributions to Higgs boson masses in the next-to-minimal supersym¬ 
metric standard model,” Phys.Lett. , vol. B314, pp. 56-63, 1993, hep-ph/9305282 

[31] T. Elliott, S. King, and P. White, “Radiative corrections to Higgs boson masses in the next-to-minimal supersym¬ 
metric Standard Model,” Phys.Rev., vol. D49, pp. 2435-2456, 1994, hep-ph/9308309 

[32] G. Yeghian, “Upper bound on the lightest Higgs mass in supersymmetric theories,” Acta Phys.Slov., vol. 49, p. 823, 
199, hep-ph/9904488 

[33] U. Ellwanger and C. Hugonie, “Masses and couplings of the lightest Higgs bosons in the (M+l) SSM,” Eur.Phys.J., 
vol. C25, pp. 297-305, 2002, hep-ph/9909260 

[34] U. Ellwanger and C. Hugonie, “Yukawa induced radiative corrections to the lightest Higgs boson mass in the 
NMSSM,” Phys.Lett., vol. B623, pp. 93-103, 2005, hep-ph/0504269 

[35] S. P. Martin, “Two loop effective potential for a general renormalizable theory and softly broken supersymmetry,” 
Phys.Rev., vol. D65, p. 116003, 2002, hep-ph/0111209 

[36] G. Degrassi and P. Slavich, “On the radiative corrections to the neutral Higgs boson masses in the NMSSM,” 
Nucl.Phys. , vol. B825, pp. 119-150, 2010, 0907.4682. 

[37] F. Staub, W. Porod, and B. Herrmann, “The Electroweak sector of the NMSSM at the one-loop level,” JHEP , 
vol. 1010, p. 040, 2010, 1007.4049. 


26 




































[38] K. Ender, T. Graf, M. Muhlleitner, and H. Rzehak, “Analysis of the NMSSM Higgs Boson Masses at One-Loop 
Level,” Phys.Rev. , vol. D85, p. 075024, 2012, 1111.4952. 

[39] T. Graf, R. Grober, M. Muhlleitner, H. Rzehak, and K. Walz, “Higgs Boson Masses in the Complex NMSSM at 
One-Loop Level,” JHEP , vol. 1210, p. 122, 2012, 1206.6806. 

[40] M. D. Goodsell, K. Nickel, and F. Staub, “Two-loop corrections to the Higgs masses in the NMSSM,” Phys.Rev. , 
vol. D91, no. 3, p. 035021, 2015, 1411.4665. 

[41] M. Muhlleitner, D. T. Nhung, H. Rzehak, and K. Walz, “Two-Loop Contributions of the Order O(atcts) to the 
Masses of the Higgs Bosons in the CP-Violating NMSSM,” JHEP , vol. 1505, p. 128, 2015, 1412.0918. 

[42] U. Ellwanger, J. F. Gunion, and C. Hugonie, “NMHDECAY: A Fortran code for the Higgs masses, couplings and 
decay widths in the NMSSM,” JHEP, vol. 0502, p. 066, 2005, hep-ph/0406215 

[43] U. Ellwanger and C. Hugonie, “NMHDECAY 2.0: An Updated program for sparticle masses, Higgs masses, 
couplings and decay widths in the NMSSM,” Comput.Phys.Commun., vol. 175, pp. 290-303, 2006, hep-ph/0508022 

[44] U. Ellwanger and C. Hugonie, “NMSPEC: A Fortran code for the sparticle and Higgs masses in the NMSSM with 
GUT scale boundary conditions,” Comput.Phys.Commun. , vol. 177, pp. 399-407, 2007, hep-ph/0612134 

[45] B. Allanach, “SOFTSUSY: a program for calculating supersymmetric spectra,” Comput.Phys.Commun., vol. 143, 
pp. 305-331, 2002, hep-ph/0104145 

[46] B. Allanach, P. Athron, L. C. Tunstall, A. Voigt, and A. Williams, “Next-to-Minimal SOFTSUSY,” 
Comput.Phys.Commun., vol. 185, pp. 2322-2339, 2014, 1311.7659. 

[47] B. C. Allanach, A. Bednyakov, and R. Ruiz de Austri, “Higher order corrections and unification in the mini¬ 
mal supersymmetric standard model: SOFTSUSY3.5,” Comput. Phys. Commun., vol. 189, pp. 192-206, 2015, 
1407.6130. 

[48] J. Baglio, R. Grober, M. Muhlleitner, D. Nhung, H. Rzehak, et al. , “NMSSMCALC: A Program Package 
for the Calculation of Loop-Corrected Higgs Boson Masses and Decay Widths in the (Complex) NMSSM,” 
Comput.Phys.Commun., vol. 185, no. 12, pp. 3372-3391, 2014, 1312.4788. 

[49] P. Athron, J.-h. Park, D. Stockinger, and A. Voigt, “FlexibleSUSY - A spectrum generator generator for super- 
symmetric models,” Comput.Phys.Commun. , vol. 190, pp. 139-172, 2015, 1406.2319. 

[50] W. Porod, “SPheno, a program for calculating supersymmetric spectra, SUSY particle decays and SUSY particle 
production at e+ e- colliders,” Comput.Phys.Commun., vol. 153, pp. 275-315, 2003, hep-ph/0301101 

[51] W. Porod and F. Staub, “SPheno 3.1: Extensions including flavour, CP-phases and models beyond the MSSM,” 

Comput.Phys.Commun. , vol. 183, pp. 2458-2469, 2012, 1104.1573. 

[52] F. Staub, “SARAH,” 2008, 0806.0538. 

[53] F. Staub, “From Superpotential to Model Files for FeynArts and CalcHep/CompHep,” Comput.Phys.Commun., 
vol. 181, pp. 1077-1086, 2010, 0909.2863. 

[54] F. Staub, “Automatic Calculation of supersymmetric Renormalization Group Equations and Self Energies,” 
Comput.Phys.Commun., vol. 182, pp. 808-833, 2011, 1002.0840. 

[55] F. Staub, “SARAH 3.2: Dirac Gauginos, UFO output, and more,” Computer Physics Communications, vol. 184, 
pp. pp. 1792-1809, 2013, 1207.0906. 

[56] F. Staub, “SARAH 4: A tool for (not only SUSY) model builders,” Comput.Phys.Commun., vol. 185, pp. 1773- 
1790, 2014, 1309.7223. 

[57] B. Allanach, A. Djouadi, J. Kneur, W. Porod, and P. Slavich, “Precise determination of the neutral Higgs boson 
masses in the MSSM,” JHEP , vol. 0409, p. 044, 2004, hep-ph/0406166 

[58] B. Allanach, S. Kraml, and W. Porod, “Theoretical uncertainties in sparticle mass predictions from computational 
tools,” JHEP , vol. 0303, p. 016, 2003, hep-ph/0302102 

[59] J. E. Kim and H. P. Nilles, “The mu Problem and the Strong CP Problem,” Phys .Lett. , vol. B138, p. 150, 1984. 

[60] S. Abel, S. Sarkar, and P. White, “On the cosmological domain wall problem for the minimally extended super- 
symmetric standard model,” Nucl.Phys. , vol. B454, pp. 663-684, 1995, hep-ph/9506359 

[61] P. Z. Skands, B. Allanach, H. Baer, C. Balazs, G. Belanger, et al., “SUSY Les Houches accord: Interfacing SUSY 
spectrum calculators, decay packages, and event generators,” JHEP, vol. 0407, p. 036, 2004, hep-ph/0311123 

[62] B. Allanach, C. Balazs, G. Belanger, M. Bernhardt, F. Boudjema, et al., “SUSY Les Houches Accord 2,” 
Comput.Phys.Commun., vol. 180, pp. 8-25, 2009, 0801.0045. 

[63] A. Djouadi, J. Kalinowski, and M. Spira, “HDECAY: A Program for Higgs boson decays in the standard model 
and its supersymmetric extension,” Comput.Phys.Commun., vol. 108, pp. 56-74, 1998, hep-ph/9704448 

[64] A. Djouadi, M. Muhlleitner, and M. Spira, “Decays of supersymmetric particles: The Program SUSY-HIT 
(SUspect-SdecaY-Hdecay-InTerface),” Acta Phys.Polon., vol. B38, pp. 635-644, 2007, hep-ph/0609292 

[65] J. M. Butterworth et_al, “THE TOOLS AND MONTE CARLO WORKING GROUP 

Summary Report from the Les Houches 2009 Workshop on TeV Colliders,” in 

Physics at TeV colliders. Proceedings, 6th Workshop, dedicated to Thomas Binoth, Les Houches, France, June 8-26, 2009, 
2010, 1003.1643. 

[66] F. Domingo, “A New Tool for the study of the CP-violating NMSSM,” JHEP , vol. 06, p. 052, 2015, 1503.07087. 

[67] A. Brignole, G. Degrassi, P. Slavich, and F. Zwirner, “On the 0(alpha(t)**2) two loop corrections to the neutral 
Higgs boson masses in the MSSM,” Nucl.Phys. , vol. B631, pp. 195-218, 2002, hep-ph/0112177 

[68] A. Dedes and P. Slavich, “Two loop corrections to radiative electroweak symmetry breaking in the MSSM,” 
Nucl.Phys., vol. B657, pp. 333-354, 2003, hep-ph/0212132 

[69] A. Dedes, G. Degrassi, and P. Slavich, “On the two loop Yukawa corrections to the MSSM Higgs boson masses at 
large tan beta,” Nucl.Phys., vol. B672, pp. 144-162, 2003, hep-ph/0305127 

[70] D. Das, U. Ellwanger, and A. M. Teixeira, “NMSDECAY: A Fortran Code for Supersymmetric Particle Decays 
in the Next-to-Minimal Supersymmetric Standard Model,” Comput.Phys.Commun. , vol. 183, pp. 774-779, 2012, 


27 







































1106.5633. 

[71] M. Muhlleitner, A. Djouadi, and Y. Mambrini, “SDECAY: A Fortran code for the decays of the supersymmetric 
particles in the MSSM,” Comput.Phys.Commun., vol. 168, pp. 46-70, 2005, hep-ph/0311167 

[72] F. Domingo and U. Ellwanger, “Updated Constraints from B Physics on the MSSM and the NMSSM,” JHEP, 
vol. 0712, p. 090, 2007, 0710.3714. 

[73] F. Domingo and U. Ellwanger, “Constraints from the Muon g-2 on the Parameter Space of the NMSSM,” JHEP , 
vol. 0807, p. 079, 2008, 0806.0733. 

[74] S. Schael et al., “Search for neutral MSSM Higgs bosons at LEP,” Eur.Phys.J., vol. C47, pp. 547-587, 2006, 
hep-ex /0602042 

[75] G. Belanger, B. Dumont, U. Ellwanger, J. Gunion, and S. Kraml, “Global fit to Higgs signal strengths and couplings 
and implications for extended Higgs sectors,” Phys.Rev., vol. D88, p. 075008, 2013, 1306.2941. 

[76] A. Delgado, G. Giudice, and P. Slavich, “Dynamical mu term in gauge mediation,” Phys.Lett., vol. B653, pp. 424- 
433, 2007, 0706.3873. 

[77] S. P. Martin and M. T. Vaughn, “Two loop renormalization group equations for soft supersymmetry breaking 
couplings,” Phys.Rev. , vol. D50, p. 2282, 1994, hep-ph/9311340 

[78] Y. Yamada, “Two loop renormalization group equations for soft SUSY breaking scalar interactions: Supergraph 
method,” Phys.Rev. , vol. D50, pp. 3537-3545, 1994, hep-ph/9401241 

[79] I. Jack and D. Jones, “Renormalization of the Fayet-Iliopoulos D term,” Phys.Lett., vol. B473, pp. 102-108, 2000, 
hep-ph/9911491 

[80] I. Jack, D. Jones, and S. Parsons, “The Fayet-Iliopoulos D term and its renormalization in softly broken super- 
symmetric theories,” Phys.Rev. , vol. D62, p. 125022, 2000, hep-ph/0007291 

[81] R. M. Fonseca, M. Malinsky, W. Porod, and F. Staub, “Running soft parameters in SUSY models with multiple 
U(l) gauge factors,” Nucl.Phys. , vol. B854, pp. 28-53, 2012, 1107.2670. 

[82] M. D. Goodsell, “Two-loop RGEs with Dirac gaugino masses,” JHEP, vol. 1301, p. 066, 2013, 1206.6697. 

[83] M. Sperling, D. Stockinger, and A. Voigt, “Renormalization of vacuum expectation values in spontaneously broken 
gauge theories,” JHEP, vol. 1307, p. 132, 2013, 1305.1548. 

[84] M. Sperling, D. Stockinger, and A. Voigt, “Renormalization of vacuum expectation values in spontaneously broken 
gauge theories: Two-loop results,” JHEP , vol. 1401, p. 068, 2014, 1310.7629. 

[85] M. D. Goodsell, K. Nickel, and F. Staub, “Two-Loop Higgs mass calculations in supersymmetric models beyond 
the MSSM with SARAH and SPheno,” Eur. Phys. J. , vol. C75, no. 1, p. 32, 2015, 1411.0675. 

[86] M. Goodsell, K. Nickel, and F. Staub, “Generic two-loop Higgs mass calculation from a diagrammatic approach,” 
Eur. Phys. J. , vol. C75, no. 6, p. 290, 2015, 1503.03098. 

[87] S. P. Martin, “Two loop scalar self energies in a general renormalizable theory at leading order in gauge couplings,” 
Phys.Rev. , vol. D70, p. 016005, 2004, hep-ph/0312092 

[88] S. P. Martin, “Evaluation of two loop selfenergy basis integrals using differential equations,” Phys.Rev., vol. D68, 
p. 075002, 2003, hep-ph/0307101 

[89] S. P. Martin, “Two-loop scalar self-energies and pole masses in a general renormalizable theory with massless gauge 
bosons,” Phys.Rev. , vol. D71, p. 116004, 2005, hep-ph/0502168 

[90] W. Porod, F. Staub, and A. Vicente, “A Flavor Kit for BSM models,” Eur. Phys. J., vol. C74, no. 8, p. 2992, 2014, 
1405.1434. 

[91] J. R. Ellis, K. Enqvist, D. V. Nanopoulos, and F. Zwirner, “Observables in Low-Energy Superstring Models,” 
Mod.Phys.Lett., vol. Al, p. 57, 1986. 

[92] R. Barbieri and G. Giudice, “Upper Bounds on Supersymmetric Particle Masses,” Nucl.Phys., vol. B306, p. 63, 
1988. 

[93] P. Bechtle, O. Brein, S. Heinemeyer, G. Weiglein, and K. E. Williams, “HiggsBounds: Confronting Arbitrary Higgs 
Sectors with Exclusion Bounds from LEP and the Tevatron,” Comput.Phys.Commun., vol. 181, pp. 138-167, 2010, 
0811.4169. 

[94] P. Bechtle, O. Brein, S. Heinemeyer, G. Weiglein, and K. E. Williams, “HiggsBounds 2.0.0: Confronting Neutral and 
Charged Higgs Sector Predictions with Exclusion Bounds from LEP and the Tevatron,” Comput.Phys.Commun., 
vol. 182, pp. 2605-2631, 2011, 1102.1898. 

[95] P. Bechtle, O. Brein, S. Heinemeyer, O. Stal, T. Stefaniak, et al., “HiggsBounds — 4: Improved Tests of Extended 
Higgs Sectors against Exclusion Bounds from LEP, the Tevatron and the LHC,” Eur.Phys.J., vol. C74, no. 3, 
p. 2693, 2014, 1311.0055. 

[96] P. Bechtle, S. Heinemeyer, O. Stal, T. Stefaniak, and G. Weiglein, “ HiggsSignals: Confronting arbitrary Higgs 
sectors with measurements at the Tevatron and the LHC,” Eur.Phys.J., vol. C74, no. 2, p. 2711, 2014, 1305.1933. 

[97] J. Alwall, M. Herquet, F. Maltoni, O. Mattelaer, and T. Stelzer, “MadGraph 5 : Going Beyond,” JHEP, vol. 1106, 
p. 128, 2011, 1106.0522. 

[98] J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, et al., “The automated computation of tree-level and 
next-to-leading order differential cross sections, and their matching to parton shower simulations,” JHEP, vol. 1407, 
p. 079, 2014, 1405.0301. 

[99] W. Kilian, T. Ohl, and J. Reuter, “WHIZARD: Simulating Multi-Particle Processes at LHC and ILC,” Eur.Phys.J., 
vol. C71, p. 1742, 2011, 0708.4233. 

[100] M. Moretti, T. Ohl, and J. Reuter, “O’Mega: An Optimizing matrix element generator,” 2001, hep-ph/0102195 

[101] F. Staub, T. Ohl, W. Porod, and C. Speckner, “A Tool Box for Implementing Supersymmetric Models,” 
Comput.Phys.Commun., vol. 183, pp. 2165-2206, 2012, 1109.5147. 

[102] V. D. Barger, M. Berger, and P. Ohmann, “The Supersymmetric particle spectrum,” Phys.Rev., vol. D49, pp. 4908- 
4930, 1994, hep-ph/9311269 


28 



































[103] S. Heinemeyer, W. Hollik, and G. Weiglein, “FeynHiggs: A Program for the calculation of the masses of the neutral 
CP even Higgs bosons in the MSSM,” Comput.Phys.Commun., vol. 124, pp. 76-89, 2000, hep-ph/9812320 

[104] S. Heinemeyer, W. Hollik, and G. Weiglein, “The Masses of the neutral CP - even Higgs bosons in the MSSM: 
Accurate analysis at the two loop level,” Eur.Phys.J., vol. C9, pp. 343-366, 1999, hep-ph/9812472 

[105] G. Degrassi, S. Heinemeyer, W. Hollik, P. Slavich, and G. Weiglein, “Towards high precision predictions for the 
MSSM Higgs sector,” Eur.Phys.J., vol. C28, pp. 133-143, 2003, hep-ph/0212020 

[106] M. Frank, T. Hahn, S. Heinemeyer, W. Hollik, H. Rzehak, et al. , “The Higgs Boson Masses and Mixings of the 
Complex MSSM in the Feynman-Diagrammatic Approach,” JHEP, vol. 0702, p. 047, 2007, hep-ph/0611326 

[107] T. Hahn, S. Heinemeyer, W. Hollik, H. Rzehak, and G. Weiglein, “High-precision predictions for the light CP-even 
Higgs Boson Mass of the MSSM,” Phys.Rev.Lett. , vol. 112, p. 141801, 2014, 1312.4937. 

[108] S. King, M. Muhlleitner, and R. Nevzorov, “NMSSM Higgs Benchmarks Near 125 GeV,” Nucl.Phys., vol. B860, 
pp. 207-244, 2012, 1201.2671. 

[109] L. J. Hall, D. Pinner, and J. T. Ruderman, “A Natural SUSY Higgs Near 126 GeV,” JHEP, vol. 1204, p. 131, 2012, 
1112.2703. 

[110] D. M. Pierce, J. A. Bagger, K. T. Matchev, and R.-j. Zhang, “Precision corrections in the minimal supersymmetric 
standard model,” Nucl.Phys. , vol. B491, pp. 3-67, 1997, hep-ph/9606211 

[111] L. Avdeev and M. Y. Kalmykov, “Pole masses of quarks in dimensional reduction,” Nucl.Phys., vol. B502, pp. 419- 
435, 1997, hep-ph/9701308 

[112] A. Bednyakov, A. Onishchenko, V. Velizhanin, and O. Veretin, “Two loop 0(alpha-s**2) MSSM corrections to the 
pole masses of heavy quarks,” Eur.Phys.J., vol. C29, pp. 87-101, 2003, hep-ph/0210258 

[113] F. Staub, “Exploring new models in all detail with SARAH,” Adv. High Energy Phys. , vol. 2015, p. 840780, 2015, 
1503.04200. 

[114] S. Fanchiotti, B. A. Kniehl, and A. Sirlin, “Incorporation of QCD effects in basic corrections of the electroweak 
theory,” Phys.Rev., vol. D48, pp. 307-331, 1993, hep-ph/9212285 

[115] G. Degrassi, S. Fanchiotti, and A. Sirlin, “Relations Between the On-shell and Ms Frameworks and the M (W) - 
M ( Z ) Interdependence,” Nucl.Phys. , vol. B351, pp. 49-69, 1991. 

[116] J. Grifols and J. Sola, “One Loop Renormalization of the Electroweak Parameters in N = 1 Supersymmetry,” 
Nucl.Phys., vol. B253, p. 47, 1985. 

[117] P. H. Chankowski, A. Dabelstein, W. Hollik, W. M. Mosle, S. Pokorski, et al., “Delta R in the MSSM,” Nucl.Phys., 
vol. B417, pp. 101-129, 1994. 

[118] C. Beskidt, W. de Boer, and D. Kazakov, “A comparison of the Higgs sectors of the CMSSM and NMSSM for a 
126 GeV Higgs boson,” Phys.Lett., vol. B726, pp. 758-766, 2013, 1308.1333. 

[119] S. P. Martin, “Three-loop Standard Model effective potential at leading order in strong and top Yukawa couplings,” 
Phys.Rev. , vol. D89, p. 013003, 2014, 1310.7553. 

[120] J. Elias-Miro, J. Espinosa, and T. Konstandin, “Taming Infrared Divergences in the Effective Potential,” JHEP, 
vol. 1408, p. 034, 2014, 1406.2652. 

[121] S. P. Martin, “Taming the Goldstone contributions to the effective potential,” Phys.Rev. , vol. D90, p. 016013, 
2014, 1406.2355. 

[122] A. Djouadi, J.-L. Kneur, and G. Moultaka, “SuSpect: A Fortran code for the supersymmetric and Higgs particle 
spectrum in the MSSM,” Comput.Phys.Commun., vol. 176, pp. 426-455, 2007, hep-ph/0211331 

[123] P. Draper, G. Lee, and C. E. M. Wagner, “Precise estimates of the Higgs mass in heavy supersymmetry,” Phys.Rev., 
vol. D89, no. 5, p. 055023, 2014, 1312.5743. 

[124] E. Bagnaschi, G. F. Giudice, P. Slavich, and A. Strumia, “Higgs Mass and Unnatural Supersymmetry,” JHEP, 
vol. 09, p. 092, 2014, 1407.4081. 

[125] J. P. Vega and G. Villadoro, “SusyHD: Higgs mass Determination in Supersymmetry,” JHEP, vol. 07, p. 159, 2015, 
1504.05200. 


29 



























